LIS Oyster Lab Tissue Processing

Updated 5/7/25 by K. Lenderman - all code up to date & graphs generated - commented out outputs for winter sampling

This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Need help? Try Stackoverflow: https://stackoverflow.com/tags/ggplot2

Attaching package: ‘dplyr’

The following objects are masked from ‘package:plyr’:

    arrange, count, desc, failwith, id, mutate, rename, summarise, summarize

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

Loading required package: lattice

Attaching package: ‘lubridate’

The following objects are masked from ‘package:base’:

    date, intersect, setdiff, union


Attaching package: ‘purrr’

The following object is masked from ‘package:plyr’:

    compact


Attaching package: ‘tidyr’

The following object is masked from ‘package:reshape2’:

    smiths

This code chunk merges all .csv files within the Tissue processing folder into one data frame and outputs the full dataset into a .csv master file. This allows us to download the raw data as a .csv, add it to the repository folder, and create the master data file without copying and pasting data in excel.

reading in .csv files from local folder

#data_all <- list.files(path = "Lab_Data_TissueProcessing/raw_data/Files_by_Month",                           # Identify all CSV files
# pattern = "*.csv", full.names = TRUE) %>% 
#lapply(read_csv) %>%   # Store all files in list
#  bind_rows          # Combine data sets into one data set 
#data_all                                            # Print data to RStudio console


#as.data.frame(data_all)  # Convert tibble to data.frame


#Filtering NAs and unnecessary columns
#data_all <- data_all %>% filter(!is.na(date_collected))
#data_all <- select(data_all, -light_regime, -oyster_zone)


#write.csv(data_all, "Master_files/tissue_processing_all_data.csv", row.names=FALSE)




##### USE THIS CODE TO MERGE DATA FILES - some of the files have columns that are not the same format (dates specifically), which was causing issues in merging. This code below should solve that problem. It converts all date columns to dates and m/d/y format #####



file_paths <- list.files(path = "Lab_Data_TissueProcessing/raw_data/Files_by_Month",
                          pattern = "*.csv", full.names = TRUE)


data_all <- lapply(file_paths, function(file_path) {
  read_csv(file_path) %>%
    mutate(
      date_collected = as.Date(date_collected,format = "%m/%d/%Y"),
      date_processed = as.Date(date_processed,format = "%m/%d/%Y"),
      date_davidsons = as.Date(date_davidsons,format = "%m/%d/%Y"),
      date_etoh = as.Date(date_etoh,format = "%m/%d/%Y")
    )
}) %>%
  bind_rows

New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 38 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, oyster_zone, condition
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): light_regime, dissection_notes, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24...
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl  (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl  (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl  (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 35 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 46 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...25, ...26, ...27...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl  (2): lab_sample, height_mm
lgl (17): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, ...23, ...24, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 998 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...25, ...26, ...27...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...25, ...26, ...27...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23,...
date  (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 36 Columns: 20── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (4): lab_id, site, condition, dissection_notes
dbl  (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (7): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (4): lab_id, site, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ......
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr   (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl   (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl  (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...26, ...27, ...28, .....
date  (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl  (2): lab_sample, height_mm
lgl (17): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...22, ...23, ...24, .....
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.

as.data.frame(data_all)  # Convert tibble to data.frame

data_all <- data_all %>% filter(!is.na(date_collected))
data_all <- select(data_all, -light_regime, -oyster_zone)


View(data_all)

write.csv(data_all, "Master_files/tissue_processing_all_data.csv", row.names=FALSE)

#A wrong value caliper input was identified for the height of sample 0923GOLD_23. This code is removing that value from the data set as we cannot conclude what this original value was. The value is 8.62. This will cause this individual to fall out of the dataset when standardized to length. This code does not completely remove the individual from the dataset.

data_all$height_mm[data_all$height_mm == "0"] <- NA

adding a month & year column to the data

data_all <- data_all %>% dplyr::mutate(date_collected= as.Date(date_collected), month = month(date_collected))

data_all <- data_all %>% dplyr::mutate(date_collected= as.Date(date_collected), year = year(date_collected))

#changing numeric month to month name
data_all$month <- factor(data_all$month, levels = c("1","2","3","4","5","6","7", "8", "9", "10", "11", "12"),
        labels=c("Jan","Feb", "March","April","May", "June", "July", "Aug", "Sept", "Oct", "Nov", "Dec"))

data_all

This chunk of code creates a numerical value in a new column for the body condition scores

data_all<- data_all %>%
 dplyr::mutate(condition_score = recode(condition, "1_very_good" = 1, "2_good" = 2, "3_good_minus"= 3, "4_fair_plus"= 4, "5_fair"= 5,"6_fair_minus"=6,"7_poor_plus"=7, "8_poor"= 8, "9_very_poor"= 9))

Warning: There was 1 warning in `dplyr::mutate()`.
ℹ In argument: `condition_score = recode(...)`.
Caused by warning:
! Unreplaced values treated as NA as `.x` is not compatible.
Please specify replacements exhaustively or supply `.default`.

head(data_all)

NA
NA

#This chunk of code is removing 0723LAUR_20 and 0723LAUR_26 from the datasheet as they have been identified as spat on shell to avoid bias in the data. During this sample collection there were animals that were significantly smaller than the single set oysters. These individuals should be removed from all monthly sampling related datasheets including disease analysis. All tissue amples will be disgarded.

data_all <- data_all %>%
  subset(lab_id != "0723LAUR_20") %>%
  subset(lab_id != "0723LAUR_26")

data_all

NA

Summary of all data - height

st_height <- summarySE(data_all%>% filter(!is.na(height_mm)), measurevar="height_mm", groupvars=c("site", "date_collected"))

st_height 

#Calculate completeness for QC
st_height$Completeness <- st_height$N /30

st_height

write.csv(st_height, "Lab_Data_TissueProcessing/output\\Completeness_tissue_processing_data.csv", row.names=FALSE)

#Mean Height

ggplot(data=data_all, aes(x=site, y=height_mm, fill=site)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Mean Shell Height ", x ="site", y = "Mean Shell Height (mm)") + facet_wrap(.~year)

#Body condition

ggplot(data=data_all, aes(x=site, y=condition_score, fill=site)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Mean Condition Score ", x ="site", y = "Mean Body Condition Score")+ scale_y_reverse()+ facet_wrap(.~year)

mean_body_condition <- data_all %>%
  dplyr::group_by(site, month, year)%>%
  dplyr::summarize(mean_bsc = mean(condition_score, na.rm = TRUE))

`summarise()` has grouped output by 'site', 'month'. You can override using the `.groups` argument.

mean_body_condition

#for overlay graph
#mean_bcs_ashc <- mean_body_condition %>%
  #filter(site == "ASHC") %>%
  #filter(year == "2024")

#Mean body condtion scores - 2023
data_2023 <- data_all%>%  filter(year=="2023")
  
ggplot(data=data_2023, aes(x=month, y=condition_score, group = month, fill = site)) +
  geom_boxplot()+ 
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = -35, vjust = 0.5, hjust=1))+
  labs(title="Mean Condition Score - 2023", x ="month", y = "Mean Body Condition Score")+ scale_y_reverse()+facet_wrap(~site)

#Mean body condtion scores - 2024
data_2024 <- data_all%>%  filter(year=="2024")
  
ggplot(data=data_2024, aes(x=month, y=condition_score, group = month, fill = site)) +
  geom_boxplot()+ 
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = -35, vjust = 0.5, hjust=1))+
  labs(title="Mean Condition Score - 2024", x ="month", y = "Mean Body Condition Score")+ scale_y_reverse()+facet_wrap(~site)

#Mean body condition scores - 2025
data_2025 <- data_all%>%  filter(year=="2025")
  
ggplot(data=data_2025, aes(x=month, y=condition_score, group = month, fill = site)) +
  geom_boxplot()+ 
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = -35, vjust = 0.5, hjust=1))+
  labs(title="Mean Condition Score - 2025", x ="month", y = "Mean Body Condition Score")+ scale_y_reverse()+facet_wrap(~site)

#Proportions graph Body condition scores

#Proportions graph Body condition scores - 2023

df_BCS_proportions2023<- data_2023 %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%

  group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%

  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.

df_BCS_proportions2023<- na.omit(df_BCS_proportions2023)

BCS_proportion_all_2023<- ggplot(data=df_BCS_proportions2023, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title="Proportion of body condition scores - 2023", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
  scale_fill_brewer(palette = "Blues", direction = -1)+
  facet_wrap(~ site)

BCS_proportion_all_2023

#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2023.pdf"), height = 7, width = 13)
#print(BCS_proportion_all_2023)
#dev.off()

#Proportions graph Body condition scores - 2024
df_BCS_proportions2024<- data_2024 %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%

  group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%

  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.

df_BCS_proportions2024<- na.omit(df_BCS_proportions2024)

BCS_proportion_all_2024<- ggplot(data=df_BCS_proportions2024, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title="Proportion of body condition scores - 2024", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
  scale_fill_brewer(palette = "Blues", direction = -1)+
  facet_wrap(~ site)

BCS_proportion_all_2024

#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2024.pdf"), height = 7, width = 13)
#print(BCS_proportion_all_2024)
#dev.off()

#Proportions graph Body condition scores - 2025
df_BCS_proportions2025<- data_2025 %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%

  group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%

  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.

df_BCS_proportions2025<- na.omit(df_BCS_proportions2025)

BCS_proportion_all_2025<- ggplot(data=df_BCS_proportions2025, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title="Proportion of body condition scores - 2025", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
  scale_fill_brewer(palette = "Blues", direction = -1)+
  facet_wrap(~ site)#, scales = "free")

BCS_proportion_all_2025

pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2025.pdf"), height = 7, width = 13)

print(BCS_proportion_all_2025)
dev.off()

png 
  2 

#Winter sampling 
# winter_bcs <- rbind(df_BCS_proportions2024, df_BCS_proportions2025)
# 
# winter_bcs <- winter_bcs %>% filter(month == "Dec"|month == "Jan"|month == "Feb")
# 
# BCS_proportion_winter<- ggplot(data=winter_bcs, aes(x=factor (month, level=c('Dec', 'Jan', 'Feb')), y= Proportion, fill=condition_score_bin)) +
#   geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+  
#    theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
#   theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
#   labs(title="Proportion of body condition scores - Winter Sampling 2024-2025", x ="Month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
#   theme(axis.text=element_text(size=12))+
#   scale_fill_brewer(palette = "Blues", direction = -1)+
#   facet_wrap(~ site)
# 
# BCS_proportion_winter

# Proportion per month of intensity scores at three sites during 2023 & 2024 & 2025
df_bcs_proportions_new <- rbind(df_BCS_proportions2023, df_BCS_proportions2024, df_BCS_proportions2025)

df_bcs_proportions_new<- df_bcs_proportions_new %>% filter(!site=="LAUR")

ggplot(data=df_bcs_proportions_new, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill",colour = "black")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(fill ="Body condition scores")+
  labs(title="Proportion of body condition scores", x ="month", y = "Proportion")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
   #scale_x_continuous("Month", breaks = c(3,4,5,6,7,8,9,10))+ 
  facet_grid(year ~site)

#Shell Pathology - all data #Includes 2023 data and all data that is in the 2024 format

shell_path <- list.files(path = "Lab_Data_TissueProcessing/raw_data/shell_pathology",                           # Identify all CSV files
 pattern = "*.csv", full.names = TRUE) %>% 
lapply(read_csv) %>%   # Store all files in list
  bind_rows          # Combine data sets into one data set 

Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (24): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discoloration, horn_add, ...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (24): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discoloration, horn_add, ...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (25): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 45 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discoloration, horn_add, ...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 15── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (14): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 28 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 30 Columns: 15── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, ...15
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolor...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (13): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 26── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (1): lab_id
lgl (25): pathology_notes, boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_dig...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 30 Columns: 14── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr  (2): lab_id, pathology_notes
lgl (12): boring_sponge, shell_scarring, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, pale_digestive, discolora...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 838 Columns: 11── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): lab_id, pathology_notes
lgl (9): boring_sponge, polydora, conchiolin_mod_severe, mud_blister, pea_crab, gill_erosion, horn_add, shell_scarring, cyst_abscess
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.

shell_path <- shell_path %>% filter(!is.na(lab_id))

shell_path <- shell_path %>% select(-c("...14":"...26"))

write.csv(shell_path, "Master_files/shell_pathology_all_data.csv", row.names=FALSE)


#Writing summary for shell pathology
shell_path_summary <- shell_path%>% separate(lab_id, c('Date_Site', 'ID'))
shell_path_summary <- select(shell_path_summary, -c("pathology_notes","ID",))
shell_path_summary$boring_sponge <- as.integer(as.logical(shell_path_summary$boring_sponge))
shell_path_summary$polydora <- as.integer(as.logical(shell_path_summary$polydora))
shell_path_summary$conchiolin_mod_severe <- as.integer(as.logical(shell_path_summary$conchiolin_mod_severe))
shell_path_summary$mud_blister <- as.integer(as.logical(shell_path_summary$mud_blister))
shell_path_summary$pea_crab <- as.integer(as.logical(shell_path_summary$pea_crab))
shell_path_summary$gill_erosion <- as.integer(as.logical(shell_path_summary$gill_erosion))
shell_path_summary$pale_digestive <- as.integer(as.logical(shell_path_summary$pale_digestive))
shell_path_summary$discoloration <- as.integer(as.logical(shell_path_summary$discoloration))
shell_path_summary$horn_add <- as.integer(as.logical(shell_path_summary$horn_add))
shell_path_summary$cyst_abscess <- as.integer(as.logical(shell_path_summary$cyst_abscess))
shell_path_summary$tumor <- as.integer(as.logical(shell_path_summary$tumor))
shell_path_summary$oyster_drill <- as.integer(as.logical(shell_path_summary$oyster_drill))
shell_path_summary$boring_snail <- as.integer(as.logical(shell_path_summary$boring_snail))
shell_path_summary$shell_scarring <- as.integer(as.logical(shell_path_summary$shell_scarring))
shell_path_count<- shell_path_summary %>%
  dplyr::group_by(Date_Site) %>%
  dplyr::summarize(Sample_count = n(),
              Boring_sponge =sum(boring_sponge),
              Polydora = sum(polydora),
              Conchiolin =sum(conchiolin_mod_severe),
              Mud_blister =sum(mud_blister),
              Pea_crab =sum(pea_crab),
              Gill_erosion =sum(gill_erosion),
              Pale_digestive =sum(pale_digestive),
              Discoloration =sum(discoloration),
              Horn =sum(horn_add),
              Cyst =sum(cyst_abscess),
              Tumor =sum(tumor),
              Oyster_drill =sum(oyster_drill),
              Boring_snail =sum(boring_snail),
              Shell_scarring =sum(shell_scarring),
              Pathogen_count = sum(boring_sponge,polydora,conchiolin_mod_severe, mud_blister, pea_crab,
                                    gill_erosion, pale_digestive, discoloration, horn_add, cyst_abscess,
                                    tumor, oyster_drill, boring_snail, shell_scarring, na.rm = TRUE)) %>%
  ungroup()
write.csv(shell_path_count, "Lab_Data_TissueProcessing/output/shell_pathology_counts.csv", row.names=FALSE)

shell_path_summary

Ash Creek Summary

df_ASHC<- data_all%>%
  filter(site=="ASHC")
df_ASHC

## Shell Height ##
st_height_ASHC <- summarySE(df_ASHC, measurevar="height_mm", groupvars=c("date_collected"))
st_height_ASHC

## Body Condition ##
#Excludes April and May due to scoring change. These months are scored categorically 'fat, medium, watery'. 

st_condition_ASHC<- summarySE(df_ASHC, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_ASHC

mean_condition_ASHC <- st_condition_ASHC %>%
  filter(!is.na(condition_score)) %>%  # Filter out rows where condition_score is NA
  summarize(mean_bcs = mean(condition_score, na.rm = TRUE))
mean_condition_ASHC

NA

ggplot(data=df_ASHC, aes(x=month, y=height_mm, group= month)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title=" Ash Creek Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)")+facet_wrap(.~year)

#Condition distribution across all sample months - 2023
df_ASHC_2023 <- df_ASHC%>%  filter(!year=="2024")

ASHC_BCS_dist2023 <-ggplot(data=df_ASHC_2023, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Ash Creek Body Condition Score - 2023", x ="Condition categorization")+
  facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2023

#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2023.pdf"), height = 7, width = 13)
#print(ASHC_BCS_dist2023)
#dev.off()

#Condition distribution across all sample months - 2024
df_ASHC_2024 <- df_ASHC%>%  filter(!year=="2023")

ASHC_BCS_dist2024 <-ggplot(data=df_ASHC_2024, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Ash Creek Body Condition Score - 2024", x ="Condition categorization")+
  facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2024

#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2024.pdf"), height = 7, width = 13)
#print(ASHC_BCS_dist2024)
#dev.off()

#Condition distribution across all sample months - 2025
df_ASHC_2025 <- df_ASHC%>%  filter(year=="2025")

ASHC_BCS_dist2025 <-ggplot(data=df_ASHC_2025, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Ash Creek Body Condition Score - 2025", x ="Condition categorization")+
  facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2025

pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2025.pdf"), height = 7, width = 13)

print(ASHC_BCS_dist2025)
dev.off()

png 
  2 

#Mean Body condition per month
ggplot(data=df_ASHC, aes(x= month, y= condition_score, group = month)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Ash Creek Mean Body Condition Score ", x ="month", y= " condition score (1-9)") + scale_y_reverse()+facet_wrap(.~year)

#ASHC Proportions graph Body condition scores - all years

ASHC_BCS_proportions<- df_ASHC %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

`summarise()` has grouped output by 'site', 'month', 'condition_score_bin'. You can override using the `.groups` argument.

ASHC_BCS_proportions<- na.omit(ASHC_BCS_proportions)

ASHC_BCS_proportions

BCS_proportion_ASHC<- ggplot(data=ASHC_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title="Ash Creek Proportion of body condition scores June- November", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
  facet_wrap(~ year)
 #scale_fill_brewer() +
  #facet_wrap(~ site)

BCS_proportion_ASHC
  
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_proportion.pdf"), height = 7, width = 13)

print(BCS_proportion_ASHC)
dev.off()

png 
  2 

% of scores >3 at ash creek - all years

df_ASHC

ASHC_precent_greater_3 <- df_ASHC %>%
  dplyr::group_by(month, site, year) %>%
  dplyr::summarise(Percentage = mean(condition_score <= 3)*100)

`summarise()` has grouped output by 'month', 'site'. You can override using the `.groups` argument.

ASHC_precent_greater_3

ASHC_BCS_percentage <- ASHC_precent_greater_3%>% 
  #filter(year =="2024")%>%
  ggplot(aes(x = month, y = Percentage)) +
  #geom_bar(width = 0.5, stat = "identity", position = "fill") +
    geom_col()+
  theme_bw() +
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
    axis.title.y = element_text(size = rel(1.3), angle = 90),
    axis.title.x = element_text(size = rel(1.3), angle = 0),
    axis.text = element_text(size = 12)
  ) +
  labs(
    title = "Ash Creek % of body condition scores >= 3",
    x = "month",
    y = "Percentage of body condition scores >= 3"
  ) +
    facet_wrap(~year)
  #scale_x_continuous( breaks = seq(5,12, by =1) )+
  #ylim(0,60)+ 
  #scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
ASHC_BCS_percentage

pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/ASHC_BCS_percentage.pdf"),height = 7, width = 13)

print(ASHC_BCS_percentage)
dev.off() 

png 
  2 

Fence Creek Summary

df_FENC<- data_all%>%
  filter(site=="FENC")
df_FENC

## Shell Height ##
st_height_FENC <- summarySE(df_FENC, measurevar="height_mm", groupvars=c("date_collected"))
st_height_FENC 

## Body Condition ##
#Excludes April and May due to scoring change. These months are scored categorically 'fat, medium, watery'. 

st_condition_FENC<- summarySE(df_FENC%>% filter(!is.na(condition_score)), measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_FENC

NA

ggplot(data=df_FENC, aes(x=month, y=height_mm, group=month)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title=" Fence Creek Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)")+ facet_wrap(.~year)

#Condition distribution across all sample months - 2023
df_FENC.2023 <- df_FENC%>%  filter(!year=="2024")

ggplot(data=df_FENC.2023, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Fence Creek Body Condition index - 2023", x ="Condition categorization")+
    facet_wrap(~ month, scales = "free")

#Condition distribution across all sample months - 2024
df_FENC.2024 <- df_FENC%>%  filter(!year=="2023")

ggplot(data=df_FENC.2024, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Fence Creek Body Condition index - 2024", x ="Condition categorization")+
    facet_wrap(~ month, scales = "free")

#Mean Body condition per month
ggplot(data=df_FENC, aes(x= month, y= condition_score, group = month)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Fence Creek Mean Body Condition Score", x ="month", y= " condition score (1-9)")+ scale_y_reverse() + facet_wrap(~ year, scales = "free")

#FENC Proportions graph Body condition scores - all years

FENC_BCS_proportions<- df_FENC %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

`summarise()` has grouped output by 'site', 'month', 'condition_score_bin'. You can override using the `.groups` argument.

FENC_BCS_proportions<- na.omit(FENC_BCS_proportions)

FENC_BCS_proportions

BCS_proportion_FENC<- ggplot(data=FENC_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title="Fence Creek Proportion of body condition scores", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
  facet_wrap(~ year)
 #scale_fill_brewer() +
  #facet_wrap(~ site)

BCS_proportion_FENC
  
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/FENC_BCS_proportion.pdf"), height = 7, width = 13)

print(BCS_proportion_FENC)
dev.off()

png 
  2 

% of scores >3 at Fence Creek - all years

#removing 0823FENC_28 bcs was NA
df_FENC <- df_FENC %>% drop_na(condition)

FENC_precent_greater_3 <- df_FENC %>%
  dplyr::group_by(month, year) %>%
  dplyr::summarise(Percentage = mean(condition_score <= 3)*100)

`summarise()` has grouped output by 'month'. You can override using the `.groups` argument.

FENC_precent_greater_3

FENC_BCS_percentage <- ggplot(FENC_precent_greater_3,aes(x = month, y = Percentage)) +
  #geom_bar(width = 0.5, stat = "identity", position = "fill") +
    geom_col()+
  theme_bw() +
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
    axis.title.y = element_text(size = rel(1.3), angle = 90),
    axis.title.x = element_text(size = rel(1.3), angle = 0),
    axis.text = element_text(size = 12)
  ) +
  labs(
    title = "Fence Creek % of body condition scores >= 3",
    x = "month",
    y = "Percentage of body condition scores >= 3"
  ) +
    facet_wrap(~year)
  #scale_x_continuous( breaks = seq(5,12, by =1) )+
  #ylim(0,60)+ 
  #scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
FENC_BCS_percentage

pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/FENC_BCS_percentage.pdf"),height = 7, width = 13)

print(FENC_BCS_percentage)
dev.off() 

png 
  2 

Gold Star Beach Summary

0524GOLD has 46 samples - we attempted to sample from 2022 and 2023 planting and consider them seperate but it appeared that due to a storm the plantings mixed together.

df_GOLD<- data_all%>%
  filter(site=="GOLD")
df_GOLD <- df_GOLD %>% filter(!is.na(height_mm))

#2023 data
df_GOLD2023 <- df_GOLD%>%  filter(!year=="2024")

#2024 data
df_GOLD2024 <- df_GOLD%>%  filter(!year=="2023")

df_GOLD$height_mm<-as.numeric(df_GOLD$height_mm)

mean_shell_height <- df_GOLD %>% mutate(year = year(date_collected)) %>% group_by(year)%>% dplyr::summarise(mean_height = mean(height_mm)) #summarySE(measurevar="height_mm", groupvars=c("year"))
mean_shell_height

## Shell Height ##
st_height_GOLD <- summarySE(df_GOLD, measurevar="height_mm", groupvars=c("date_collected"))
st_height_GOLD 

## Body Condition ##
st_condition_GOLD<- summarySE(df_GOLD, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_GOLD

NA

#GOLD Proportions graph Body condition scores

GOLD_BCS_proportions<- df_GOLD %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, year, condition_score_bin) %>% dplyr::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.

GOLD_BCS_proportions<- na.omit(GOLD_BCS_proportions)

BCS_proportion_GOLD <- ggplot(data=GOLD_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title="Gold Star Beach Proportion of body condition scores", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
  facet_wrap(~ year)
BCS_proportion_GOLD


  
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/GOLD_BCS_proportion.pdf"), height = 7, width = 13)

print(BCS_proportion_GOLD)
dev.off()

png 
  2 

% of scores >3 at gold june - november

df_GOLD

GOLD_precent_greater_3 <- df_GOLD %>%
  dplyr::group_by(month, site, year) %>%
  dplyr::summarise(Percentage = mean(condition_score <= 3)*100)

`summarise()` has grouped output by 'month', 'site'. You can override using the `.groups` argument.

GOLD_precent_greater_3

GOLD_BCS_percentage <- GOLD_precent_greater_3 %>%
  ggplot(aes(x = month, y = Percentage)) +
  #geom_bar(width = 0.5, stat = "identity", position = "fill") +
    geom_col()+
  theme_bw() +
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
    axis.title.y = element_text(size = rel(1.3), angle = 90),
    axis.title.x = element_text(size = rel(1.3), angle = 0),
    axis.text = element_text(size = 12)
  ) +
  labs(
    title = "Gold Star Beach % of body condition scores >= 3",
    x = "month",
    y = "Percentage of body condition scores >= 3"
  ) +facet_wrap(~year)
  #+ scale_x_continuous( breaks = seq(5,12, by =1) )+
  #ylim(0,60)+ 
  #scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10)) 

    
GOLD_BCS_percentage

pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/GOLD_BCS_percentage.pdf"),height = 7, width = 13)

print(GOLD_BCS_percentage)
dev.off() 

png 
  2 

ggplot(data=df_GOLD, aes(x=month, y=height_mm, group= month)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Goldstar Beach Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)") + #scale_x_continuous(limits= c(4,11), breaks = seq(5,10, by =1))+
  geom_smooth(method = "lm", se = FALSE)+ facet_wrap(~year)

#Condition distribution across all sample months - 2023
ggplot(data=df_GOLD2023, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Goldstar Beach  Body Condition Score - 2023 ", x ="Condition categorization")+
    facet_wrap(~ month, scales = "free")

#Condition distribution across all sample months - 2024
ggplot(data=df_GOLD2024, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title="Goldstar Beach  Body Condition Score - 2024", x ="Condition categorization")+
    facet_wrap(~ month, scales = "free")

#Mean Body condition per month
ggplot(data=df_GOLD, aes(x= month, y= condition_score, group = month)) +
  geom_boxplot()+  #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  labs(title="Goldstar Beach  Mean Body Condition Score ", x ="month", y= " condition score (1-9)")+ scale_y_reverse() + facet_wrap(~year)

  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+

Laurel Hollow Summary - NO LONGER A STUDY SITE

<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuZGZfTEFVUjwtIGRhdGFfYWxsJT4lXG4gIGZpbHRlcihzaXRlPT1cIkxBVVJcIilcblxuIyMgU2hlbGwgSGVpZ2h0ICMjXG5zdF9oZWlnaHRfTEFVUiA8LSBzdW1tYXJ5U0UoZGZfTEFVUiwgbWVhc3VyZXZhcj1cImhlaWdodF9tbVwiLCBncm91cHZhcnM9YyhcImRhdGVfY29sbGVjdGVkXCIpKVxuc3RfaGVpZ2h0X0xBVVIgXG5cbiMjIEJvZHkgQ29uZGl0aW9uICMjXG4jIyBFeGNsdWRlcyBNYXkgZHVlIHRvIHNjb3Jpbmcgc3lzdGVtIGNoYW5nZS4gTWF5IHNjb3JlZCBmYXQsIG1lZGl1bSwgd2F0ZXJ5LiBcbnN0X2NvbmRpdGlvbl9MQVVSPC0gc3VtbWFyeVNFKGRmX0xBVVIsIG1lYXN1cmV2YXIgPSBcImNvbmRpdGlvbl9zY29yZVwiLCBncm91cHZhcnMgPSBjKFwiZGF0ZV9jb2xsZWN0ZWRcIikpXG5zdF9jb25kaXRpb25fTEFVUlxuYGBgIn0= -->

```r
df_LAUR<- data_all%>%
  filter(site==\LAUR\)

## Shell Height ##
st_height_LAUR <- summarySE(df_LAUR, measurevar=\height_mm\, groupvars=c(\date_collected\))
st_height_LAUR 

## Body Condition ##
## Excludes May due to scoring system change. May scored fat, medium, watery. 
st_condition_LAUR<- summarySE(df_LAUR, measurevar = \condition_score\, groupvars = c(\date_collected\))
st_condition_LAUR
```

<!-- rnb-source-end -->

```r
df_LAUR<- data_all%>%
  filter(site==\LAUR\)

## Shell Height ##
st_height_LAUR <- summarySE(df_LAUR, measurevar=\height_mm\, groupvars=c(\date_collected\))
st_height_LAUR 

## Body Condition ##
## Excludes May due to scoring system change. May scored fat, medium, watery. 
st_condition_LAUR<- summarySE(df_LAUR, measurevar = \condition_score\, groupvars = c(\date_collected\))
st_condition_LAUR

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-text-begin -->



<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVaMmR3Ykc5MEtHUmhkR0U5WkdaZlRFRlZVaXdnWVdWektIZzliVzl1ZEdnc0lIazlhR1ZwWjJoMFgyMXRMQ0JuY205MWNEMGdiVzl1ZEdncEtTQXJYRzRnSUdkbGIyMWZZbTk0Y0d4dmRDZ3BLeUFnSTNOallXeGxYMlpwYkd4ZmJXRnVkV0ZzS0haaGJIVmxjejFqS0Z3aVptOXlaWE4wWjNKbFpXNWNJaXhjSW05eVlXNW5aVndpTENCY0luQjFjbkJzWlZ3aUtTa3JYRzRnSUhSb1pXMWxYMkozS0NrZ0t5QWdkR2hsYldVb2NHRnVaV3d1WjNKcFpDNXRZV3B2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1N3Z2NHRnVaV3d1WjNKcFpDNXRhVzV2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1NrcklGeHVJQ0FqZEdobGJXVW9ZWGhwY3k1MFpYaDBMbmdnUFNCbGJHVnRaVzUwWDNSbGVIUW9ZVzVuYkdVZ1BTQTVNQ3dnZG1wMWMzUWdQU0F3TGpVc0lHaHFkWE4wUFRFcEtTdGNiaUFnYkdGaWN5aDBhWFJzWlQxY0lpQk1ZWFZ5Wld3Z1NHOXNiRzkzSUUxbFlXNGdVMmhsYkd3Z1NHVnBaMmgwSUMwZ01qQXlNMXdpTENCNElEMWNJbTF2Ym5Sb1hDSXNJSGtnUFNCY0lrMWxZVzRnVTJobGJHd2dTR1ZwWjJoMElDaHRiU2xjSWlsY2JseHVYRzVnWUdBaWZRPT0gLS0+XG5cbmBgYHJcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PW1vbnRoLCB5PWhlaWdodF9tbSwgZ3JvdXA9IG1vbnRoKSkgK1xuICBnZW9tX2JveHBsb3QoKSsgICNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YyhcXGZvcmVzdGdyZWVuXFwsXFxvcmFuZ2VcXCwgXFxwdXJwbGVcXCkpK1xuICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgI3RoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gMC41LCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFwgTGF1cmVsIEhvbGxvdyBNZWFuIFNoZWxsIEhlaWdodCAtIDIwMjNcXCwgeCA9XFxtb250aFxcLCB5ID0gXFxNZWFuIFNoZWxsIEhlaWdodCAobW0pXFwpXG5cblxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->

ggplot(data=df_LAUR, aes(x=month, y=height_mm, group= month)) +
  geom_boxplot()+  #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title=\ Laurel Hollow Mean Shell Height - 2023\, x =\month\, y = \Mean Shell Height (mm)\)

<!-- rnb-output-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dVoyZHdiRzkwS0dSaGRHRTlaR1pmVEVGVlVpd2dZV1Z6S0hnOWJXOXVkR2dzSUhrOWFHVnBaMmgwWDIxdExDQm5jbTkxY0QwZ2JXOXVkR2dwS1NBclhHNGdJR2RsYjIxZlltOTRjR3h2ZENncEt5QWdJM05qWVd4bFgyWnBiR3hmYldGdWRXRnNLSFpoYkhWbGN6MWpLRnhjWm05eVpYTjBaM0psWlc1Y1hDeGNYRzl5WVc1blpWeGNMQ0JjWEhCMWNuQnNaVnhjS1NrclhHNGdJSFJvWlcxbFgySjNLQ2tnS3lBZ2RHaGxiV1VvY0dGdVpXd3VaM0pwWkM1dFlXcHZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU3dnY0dGdVpXd3VaM0pwWkM1dGFXNXZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU2tySUZ4dUlDQWpkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBNU1Dd2dkbXAxYzNRZ1BTQXdMalVzSUdocWRYTjBQVEVwS1N0Y2JpQWdiR0ZpY3loMGFYUnNaVDFjWENCTVlYVnlaV3dnU0c5c2JHOTNJRTFsWVc0Z1UyaGxiR3dnU0dWcFoyaDBJQzBnTWpBeU0xeGNMQ0I0SUQxY1hHMXZiblJvWEZ3c0lIa2dQU0JjWEUxbFlXNGdVMmhsYkd3Z1NHVnBaMmgwSUNodGJTbGNYQ2xjYmx4dVhHNWdZR0JjYm1CZ1lDSjkgLS0+XG5cbmBgYHJcbmBgYHJcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PW1vbnRoLCB5PWhlaWdodF9tbSwgZ3JvdXA9IG1vbnRoKSkgK1xuICBnZW9tX2JveHBsb3QoKSsgICNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YyhcXGZvcmVzdGdyZWVuXFwsXFxvcmFuZ2VcXCwgXFxwdXJwbGVcXCkpK1xuICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgI3RoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gMC41LCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFwgTGF1cmVsIEhvbGxvdyBNZWFuIFNoZWxsIEhlaWdodCAtIDIwMjNcXCwgeCA9XFxtb250aFxcLCB5ID0gXFxNZWFuIFNoZWxsIEhlaWdodCAobW0pXFwpXG5cblxuYGBgXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuZ2dwbG90KGRhdGE9ZGZfTEFVUiwgYWVzKHg9bW9udGgsIHk9aGVpZ2h0X21tLCBncm91cD0gbW9udGgpKSArXG4gIGdlb21fYm94cGxvdCgpKyAgI3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcz1jKFxcZm9yZXN0Z3JlZW5cXCxcXG9yYW5nZVxcLCBcXHB1cnBsZVxcKSkrXG4gIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbiAgbGFicyh0aXRsZT1cXCBMYXVyZWwgSG9sbG93IE1lYW4gU2hlbGwgSGVpZ2h0IC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXE1lYW4gU2hlbGwgSGVpZ2h0IChtbSlcXClcblxuXG5gYGBcbmBgYCJ9 -->

```r
```r
ggplot(data=df_LAUR, aes(x=month, y=height_mm, group= month)) +
  geom_boxplot()+  #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title=\ Laurel Hollow Mean Shell Height - 2023\, x =\month\, y = \Mean Shell Height (mm)\)

```
```

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVJME52Ym1ScGRHbHZiaUJrYVhOMGNtbGlkWFJwYjI0Z1lXTnliM056SUdGc2JDQnpZVzF3YkdVZ2JXOXVkR2h6WEc1blozQnNiM1FvWkdGMFlUMWtabDlNUVZWU0xDQmhaWE1vZUQwZ1kyOXVaR2wwYVc5dUtTa2dLMXh1SUNCblpXOXRYMkpoY2lncEt5QWdJM05qWVd4bFgyWnBiR3hmYldGdWRXRnNLSFpoYkhWbGN6MWpLRndpWm05eVpYTjBaM0psWlc1Y0lpeGNJbTl5WVc1blpWd2lMQ0JjSW5CMWNuQnNaVndpS1NrclhHNGdJSFJvWlcxbFgySjNLQ2tnS3lBZ2RHaGxiV1VvY0dGdVpXd3VaM0pwWkM1dFlXcHZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU3dnY0dGdVpXd3VaM0pwWkM1dGFXNXZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU2tySUZ4dUlDQWpkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBNU1Dd2dkbXAxYzNRZ1BTQXdMalVzSUdocWRYTjBQVEVwS1N0Y2JpQWdiR0ZpY3loMGFYUnNaVDFjSWt4aGRYSmxiQ0JJYjJ4c2IzY2dRbTlrZVNCRGIyNWthWFJwYjI0Z2FXNWtaWGdnTFNBeU1ESXpYQ0lzSUhnZ1BWd2lRMjl1WkdsMGFXOXVJR05oZEdWbmIzSnBlbUYwYVc5dVhDSXBLMXh1SUNBZ0lHWmhZMlYwWDNkeVlYQW9maUJ0YjI1MGFDd2djMk5oYkdWeklEMGdYQ0ptY21WbFhDSXBYRzVjYmx4dUkwMWxZVzRnUW05a2VTQmpiMjVrYVhScGIyNGdjR1Z5SUcxdmJuUm9YRzVuWjNCc2IzUW9aR0YwWVQxa1psOU1RVlZTTENCaFpYTW9lRDBnYlc5dWRHZ3NJSGs5SUdOdmJtUnBkR2x2Ymw5elkyOXlaU3dnWjNKdmRYQWdQU0J0YjI1MGFDa3BJQ3RjYmlBZ1oyVnZiVjlpYjNod2JHOTBLQ2tySUNBamMyTmhiR1ZmWm1sc2JGOXRZVzUxWVd3b2RtRnNkV1Z6UFdNb1hDSm1iM0psYzNSbmNtVmxibHdpTEZ3aWIzSmhibWRsWENJc0lGd2ljSFZ5Y0d4bFhDSXBLU3RjYmlBZ2RHaGxiV1ZmWW5jb0tTQXJJQ0IwYUdWdFpTaHdZVzVsYkM1bmNtbGtMbTFoYW05eUlEMGdaV3hsYldWdWRGOWliR0Z1YXlncExDQndZVzVsYkM1bmNtbGtMbTFwYm05eUlEMGdaV3hsYldWdWRGOWliR0Z1YXlncEtTc2dYRzRnSUd4aFluTW9kR2wwYkdVOVhDSk1ZWFZ5Wld3Z1NHOXNiRzkzSUUxbFlXNGdRbTlrZVNCRGIyNWthWFJwYjI0Z1UyTnZjbVVnTFNBeU1ESXpYQ0lzSUhnZ1BWd2liVzl1ZEdoY0lpd2dlVDBnWENJZ1kyOXVaR2wwYVc5dUlITmpiM0psSUNneExUa3BYQ0lwS3lCelkyRnNaVjk1WDNKbGRtVnljMlVvS1Z4dUlDQWpkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBNU1Dd2dkbXAxYzNRZ1BTQXdMalVzSUdocWRYTjBQVEVwS1N0Y2JtQmdZQ0o5IC0tPlxuXG5gYGByXG4jQ29uZGl0aW9uIGRpc3RyaWJ1dGlvbiBhY3Jvc3MgYWxsIHNhbXBsZSBtb250aHNcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PSBjb25kaXRpb24pKSArXG4gIGdlb21fYmFyKCkrICAjc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoXFxmb3Jlc3RncmVlblxcLFxcb3JhbmdlXFwsIFxccHVycGxlXFwpKStcbiAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gICN0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCB2anVzdCA9IDAuNSwgaGp1c3Q9MSkpK1xuICBsYWJzKHRpdGxlPVxcTGF1cmVsIEhvbGxvdyBCb2R5IENvbmRpdGlvbiBpbmRleCAtIDIwMjNcXCwgeCA9XFxDb25kaXRpb24gY2F0ZWdvcml6YXRpb25cXCkrXG4gICAgZmFjZXRfd3JhcCh+IG1vbnRoLCBzY2FsZXMgPSBcXGZyZWVcXClcblxuXG4jTWVhbiBCb2R5IGNvbmRpdGlvbiBwZXIgbW9udGhcbmdncGxvdChkYXRhPWRmX0xBVVIsIGFlcyh4PSBtb250aCwgeT0gY29uZGl0aW9uX3Njb3JlLCBncm91cCA9IG1vbnRoKSkgK1xuICBnZW9tX2JveHBsb3QoKSsgICNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9YyhcXGZvcmVzdGdyZWVuXFwsXFxvcmFuZ2VcXCwgXFxwdXJwbGVcXCkpK1xuICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgbGFicyh0aXRsZT1cXExhdXJlbCBIb2xsb3cgTWVhbiBCb2R5IENvbmRpdGlvbiBTY29yZSAtIDIwMjNcXCwgeCA9XFxtb250aFxcLCB5PSBcXCBjb25kaXRpb24gc2NvcmUgKDEtOSlcXCkrIHNjYWxlX3lfcmV2ZXJzZSgpXG4gICN0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCB2anVzdCA9IDAuNSwgaGp1c3Q9MSkpK1xuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->

#Condition distribution across all sample months
ggplot(data=df_LAUR, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title=\Laurel Hollow Body Condition index - 2023\, x =\Condition categorization\)+
    facet_wrap(~ month, scales = \free\)


#Mean Body condition per month
ggplot(data=df_LAUR, aes(x= month, y= condition_score, group = month)) +
  geom_boxplot()+  #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  labs(title=\Laurel Hollow Mean Body Condition Score - 2023\, x =\month\, y= \ condition score (1-9)\)+ scale_y_reverse()
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+

<!-- rnb-output-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dUkwTnZibVJwZEdsdmJpQmthWE4wY21saWRYUnBiMjRnWVdOeWIzTnpJR0ZzYkNCellXMXdiR1VnYlc5dWRHaHpYRzVuWjNCc2IzUW9aR0YwWVQxa1psOU1RVlZTTENCaFpYTW9lRDBnWTI5dVpHbDBhVzl1S1NrZ0sxeHVJQ0JuWlc5dFgySmhjaWdwS3lBZ0kzTmpZV3hsWDJacGJHeGZiV0Z1ZFdGc0tIWmhiSFZsY3oxaktGeGNabTl5WlhOMFozSmxaVzVjWEN4Y1hHOXlZVzVuWlZ4Y0xDQmNYSEIxY25Cc1pWeGNLU2tyWEc0Z0lIUm9aVzFsWDJKM0tDa2dLeUFnZEdobGJXVW9jR0Z1Wld3dVozSnBaQzV0WVdwdmNpQTlJR1ZzWlcxbGJuUmZZbXhoYm1zb0tTd2djR0Z1Wld3dVozSnBaQzV0YVc1dmNpQTlJR1ZzWlcxbGJuUmZZbXhoYm1zb0tTa3JJRnh1SUNBamRHaGxiV1VvWVhocGN5NTBaWGgwTG5nZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvWVc1bmJHVWdQU0E1TUN3Z2RtcDFjM1FnUFNBd0xqVXNJR2hxZFhOMFBURXBLU3RjYmlBZ2JHRmljeWgwYVhSc1pUMWNYRXhoZFhKbGJDQkliMnhzYjNjZ1FtOWtlU0JEYjI1a2FYUnBiMjRnYVc1a1pYZ2dMU0F5TURJelhGd3NJSGdnUFZ4Y1EyOXVaR2wwYVc5dUlHTmhkR1ZuYjNKcGVtRjBhVzl1WEZ3cEsxeHVJQ0FnSUdaaFkyVjBYM2R5WVhBb2ZpQnRiMjUwYUN3Z2MyTmhiR1Z6SUQwZ1hGeG1jbVZsWEZ3cFhHNWNibHh1STAxbFlXNGdRbTlrZVNCamIyNWthWFJwYjI0Z2NHVnlJRzF2Ym5Sb1hHNW5aM0JzYjNRb1pHRjBZVDFrWmw5TVFWVlNMQ0JoWlhNb2VEMGdiVzl1ZEdnc0lIazlJR052Ym1ScGRHbHZibDl6WTI5eVpTd2daM0p2ZFhBZ1BTQnRiMjUwYUNrcElDdGNiaUFnWjJWdmJWOWliM2h3Ykc5MEtDa3JJQ0FqYzJOaGJHVmZabWxzYkY5dFlXNTFZV3dvZG1Gc2RXVnpQV01vWEZ4bWIzSmxjM1JuY21WbGJseGNMRnhjYjNKaGJtZGxYRndzSUZ4Y2NIVnljR3hsWEZ3cEtTdGNiaUFnZEdobGJXVmZZbmNvS1NBcklDQjBhR1Z0WlNod1lXNWxiQzVuY21sa0xtMWhhbTl5SUQwZ1pXeGxiV1Z1ZEY5aWJHRnVheWdwTENCd1lXNWxiQzVuY21sa0xtMXBibTl5SUQwZ1pXeGxiV1Z1ZEY5aWJHRnVheWdwS1NzZ1hHNGdJR3hoWW5Nb2RHbDBiR1U5WEZ4TVlYVnlaV3dnU0c5c2JHOTNJRTFsWVc0Z1FtOWtlU0JEYjI1a2FYUnBiMjRnVTJOdmNtVWdMU0F5TURJelhGd3NJSGdnUFZ4Y2JXOXVkR2hjWEN3Z2VUMGdYRndnWTI5dVpHbDBhVzl1SUhOamIzSmxJQ2d4TFRrcFhGd3BLeUJ6WTJGc1pWOTVYM0psZG1WeWMyVW9LVnh1SUNBamRHaGxiV1VvWVhocGN5NTBaWGgwTG5nZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvWVc1bmJHVWdQU0E1TUN3Z2RtcDFjM1FnUFNBd0xqVXNJR2hxZFhOMFBURXBLU3RjYm1CZ1lGeHVZR0JnSW4wPSAtLT5cblxuYGBgclxuYGBgclxuI0NvbmRpdGlvbiBkaXN0cmlidXRpb24gYWNyb3NzIGFsbCBzYW1wbGUgbW9udGhzXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gY29uZGl0aW9uKSkgK1xuICBnZW9tX2JhcigpKyAgI3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcz1jKFxcZm9yZXN0Z3JlZW5cXCxcXG9yYW5nZVxcLCBcXHB1cnBsZVxcKSkrXG4gIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbiAgbGFicyh0aXRsZT1cXExhdXJlbCBIb2xsb3cgQm9keSBDb25kaXRpb24gaW5kZXggLSAyMDIzXFwsIHggPVxcQ29uZGl0aW9uIGNhdGVnb3JpemF0aW9uXFwpK1xuICAgIGZhY2V0X3dyYXAofiBtb250aCwgc2NhbGVzID0gXFxmcmVlXFwpXG5cblxuI01lYW4gQm9keSBjb25kaXRpb24gcGVyIG1vbnRoXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gbW9udGgsIHk9IGNvbmRpdGlvbl9zY29yZSwgZ3JvdXAgPSBtb250aCkpICtcbiAgZ2VvbV9ib3hwbG90KCkrICAjc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoXFxmb3Jlc3RncmVlblxcLFxcb3JhbmdlXFwsIFxccHVycGxlXFwpKStcbiAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IE1lYW4gQm9keSBDb25kaXRpb24gU2NvcmUgLSAyMDIzXFwsIHggPVxcbW9udGhcXCwgeT0gXFwgY29uZGl0aW9uIHNjb3JlICgxLTkpXFwpKyBzY2FsZV95X3JldmVyc2UoKVxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbmBgYFxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuI0NvbmRpdGlvbiBkaXN0cmlidXRpb24gYWNyb3NzIGFsbCBzYW1wbGUgbW9udGhzXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gY29uZGl0aW9uKSkgK1xuICBnZW9tX2JhcigpKyAgI3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcz1jKFxcZm9yZXN0Z3JlZW5cXCxcXG9yYW5nZVxcLCBcXHB1cnBsZVxcKSkrXG4gIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbiAgbGFicyh0aXRsZT1cXExhdXJlbCBIb2xsb3cgQm9keSBDb25kaXRpb24gaW5kZXggLSAyMDIzXFwsIHggPVxcQ29uZGl0aW9uIGNhdGVnb3JpemF0aW9uXFwpK1xuICAgIGZhY2V0X3dyYXAofiBtb250aCwgc2NhbGVzID0gXFxmcmVlXFwpXG5cblxuI01lYW4gQm9keSBjb25kaXRpb24gcGVyIG1vbnRoXG5nZ3Bsb3QoZGF0YT1kZl9MQVVSLCBhZXMoeD0gbW9udGgsIHk9IGNvbmRpdGlvbl9zY29yZSwgZ3JvdXAgPSBtb250aCkpICtcbiAgZ2VvbV9ib3hwbG90KCkrICAjc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoXFxmb3Jlc3RncmVlblxcLFxcb3JhbmdlXFwsIFxccHVycGxlXFwpKStcbiAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IE1lYW4gQm9keSBDb25kaXRpb24gU2NvcmUgLSAyMDIzXFwsIHggPVxcbW9udGhcXCwgeT0gXFwgY29uZGl0aW9uIHNjb3JlICgxLTkpXFwpKyBzY2FsZV95X3JldmVyc2UoKVxuICAjdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAwLjUsIGhqdXN0PTEpKStcbmBgYFxuYGBgIn0= -->

```r
```r
#Condition distribution across all sample months
ggplot(data=df_LAUR, aes(x= condition)) +
  geom_bar()+  #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
  labs(title=\Laurel Hollow Body Condition index - 2023\, x =\Condition categorization\)+
    facet_wrap(~ month, scales = \free\)


#Mean Body condition per month
ggplot(data=df_LAUR, aes(x= month, y= condition_score, group = month)) +
  geom_boxplot()+  #scale_fill_manual(values=c(\forestgreen\,\orange\, \purple\))+
  theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  labs(title=\Laurel Hollow Mean Body Condition Score - 2023\, x =\month\, y= \ condition score (1-9)\)+ scale_y_reverse()
  #theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
```
```

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-text-begin -->

#LAUR Proportions graph Body condition scores - 2023 (only sampled one year)


<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVYRzVNUVZWU1gwSkRVMTl3Y205d2IzSjBhVzl1Y3p3dElHUm1YMHhCVlZJZ0pUNGxYRzRnSUcxMWRHRjBaU2hqYjI1a2FYUnBiMjVmYzJOdmNtVmZiblZ0WlhKcFl5QTlJR0Z6TG01MWJXVnlhV01vWTI5dVpHbDBhVzl1WDNOamIzSmxLU3hqYjI1a2FYUnBiMjVmYzJOdmNtVmZZbWx1SUQwZ1kyRnpaVjkzYUdWdUtHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQXhJSDRnWENJeFhDSXNJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0F5SUg1Y0lqSmNJaXhqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTXlCK0lGd2lNMXdpTENCamIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ05DQitYQ0kwWENJc1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlEVWdmaUJjSWpWY0lpd2dZMjl1WkdsMGFXOXVYM05qYjNKbElEMDlJRFlnZmx3aU5sd2lMR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0EzSUg0Z1hDSTNYQ0lzSUdOdmJtUnBkR2x2Ymw5elkyOXlaU0E5UFNBNElINWNJamhjSWl3Z1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOU9TQitYQ0k1WENJc0lGUlNWVVVnZmlCaGN5NWphR0Z5WVdOMFpYSW9ZMjl1WkdsMGFXOXVYM05qYjNKbEtTa3BJQ1UrSlZ4dUlDQm5jbTkxY0Y5aWVTaHphWFJsTENCdGIyNTBhQ3dnWTI5dVpHbDBhVzl1WDNOamIzSmxYMkpwYml3Z2VXVmhjaWtnSlQ0bElHUndiSGx5T2pwemRXMXRZWEpwYzJVb1EyOTFiblE5SUc0b0tTa2dKVDRsWEc0Z0lIVnVaM0p2ZFhBb0tTQWxQaVZjYmlBZ2JYVjBZWFJsS0ZCeWIzQnZjblJwYjI0Z1BTQkRiM1Z1ZEM5emRXMG9RMjkxYm5RcEtWeHVYRzVNUVZWU1gwSkRVMTl3Y205d2IzSjBhVzl1Y3p3dElHNWhMbTl0YVhRb1RFRlZVbDlDUTFOZmNISnZjRzl5ZEdsdmJuTXBYRzVjYmt4QlZWSmZRa05UWDNCeWIzQnZjblJwYjI1elhHNWNia0pEVTE5d2NtOXdiM0owYVc5dVgweEJWVkk4TFNCblozQnNiM1FvWkdGMFlUMU1RVlZTWDBKRFUxOXdjbTl3YjNKMGFXOXVjeXdnWVdWektIZzliVzl1ZEdnc0lIazlJRkJ5YjNCdmNuUnBiMjRzSUdacGJHdzlZMjl1WkdsMGFXOXVYM05qYjNKbFgySnBiaWtwSUN0Y2JpQWdaMlZ2YlY5aVlYSW9kMmxrZEdnZ1BTQXVOU3dnYzNSaGREMWNJbWxrWlc1MGFYUjVYQ0lzSUhCdmMybDBhVzl1SUQwZ1hDSm1hV3hzWENJcEt5QWdYRzRnSUNCMGFHVnRaVjlpZHlncElDc2dJSFJvWlcxbEtIQmhibVZzTG1keWFXUXViV0ZxYjNJZ1BTQmxiR1Z0Wlc1MFgySnNZVzVyS0Nrc0lIQmhibVZzTG1keWFXUXViV2x1YjNJZ1BTQmxiR1Z0Wlc1MFgySnNZVzVyS0NrcEt5QmNiaUFnZEdobGJXVW9ZWGhwY3k1MFpYaDBMbmdnUFNCbGJHVnRaVzUwWDNSbGVIUW9ZVzVuYkdVZ1BTQTBOU3dnZG1wMWMzUWdQU0F4TENCb2FuVnpkRDB4S1NrclhHNGdJR3hoWW5Nb2RHbDBiR1U5WENKTVlYVnlaV3dnU0c5c2JHOTNJRkJ5YjNCdmNuUnBiMjRnYjJZZ1ltOWtlU0JqYjI1a2FYUnBiMjRnYzJOdmNtVnpJQzBnTWpBeU0xd2lMQ0I0SUQxY0ltMXZiblJvWENJc0lIa2dQU0JjSWxCeWIzQnZjblJwYjI0Z2IyWWdZbTlrZVNCamIyNWthWFJwYjI0Z2MyTnZjbVZ6SUZ3aUtTc2dkR2hsYldVb1lYaHBjeTUwYVhSc1pTNTVJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFRrd0tTd2dZWGhwY3k1MGFYUnNaUzU0SUQwZ1pXeGxiV1Z1ZEY5MFpYaDBLSE5wZW1VZ1BTQnlaV3dvTVM0ektTd2dZVzVuYkdVZ1BTQXdLU2tyWEc0Z0lIUm9aVzFsS0dGNGFYTXVkR1Y0ZEQxbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlQweE1pa3BYRzRnSUNObVlXTmxkRjkzY21Gd0tINGdlV1ZoY2lsY2JpQWpjMk5oYkdWZlptbHNiRjlpY21WM1pYSW9LU0FyWEc0Z0lDTm1ZV05sZEY5M2NtRndLSDRnYzJsMFpTbGNibHh1UWtOVFgzQnliM0J2Y25ScGIyNWZURUZWVWx4dUlDQmNibkJrWmlod1lYTjBaVEFvY0dGMGFDQTlJRndpVEdGaVgwUmhkR0ZmVkdsemMzVmxVSEp2WTJWemMybHVaeTl2ZFhSd2RYUmNJaUFzWENJdlRFRlZVbDlDUTFOZmNISnZjRzl5ZEdsdmJpNXdaR1pjSWlrc0lHaGxhV2RvZENBOUlEY3NJSGRwWkhSb0lEMGdNVE1wWEc1d2NtbHVkQ2hDUTFOZmNISnZjRzl5ZEdsdmJsOU1RVlZTS1Z4dVpHVjJMbTltWmlncFhHNWdZR0FpZlE9PSAtLT5cblxuYGBgclxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIGRmX0xBVVIgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbiwgeWVhcikgJT4lIGRwbHlyOjpzdW1tYXJpc2UoQ291bnQ9IG4oKSkgJT4lXG4gIHVuZ3JvdXAoKSAlPiVcbiAgbXV0YXRlKFByb3BvcnRpb24gPSBDb3VudC9zdW0oQ291bnQpKVxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoTEFVUl9CQ1NfcHJvcG9ydGlvbnMpXG5cbkxBVVJfQkNTX3Byb3BvcnRpb25zXG5cbkJDU19wcm9wb3J0aW9uX0xBVVI8LSBnZ3Bsb3QoZGF0YT1MQVVSX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpKyAgXG4gICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgdmp1c3QgPSAxLCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIFxcKSsgdGhlbWUoYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPTkwKSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSkrXG4gIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMikpXG4gICNmYWNldF93cmFwKH4geWVhcilcbiAjc2NhbGVfZmlsbF9icmV3ZXIoKSArXG4gICNmYWNldF93cmFwKH4gc2l0ZSlcblxuQkNTX3Byb3BvcnRpb25fTEFVUlxuICBcbnBkZihwYXN0ZTAocGF0aCA9IFxcTGFiX0RhdGFfVGlzc3VlUHJvY2Vzc2luZy9vdXRwdXRcXCAsXFwvTEFVUl9CQ1NfcHJvcG9ydGlvbi5wZGZcXCksIGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChCQ1NfcHJvcG9ydGlvbl9MQVVSKVxuZGV2Lm9mZigpXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->

LAUR_BCS_proportions<- df_LAUR %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

LAUR_BCS_proportions<- na.omit(LAUR_BCS_proportions)

LAUR_BCS_proportions

BCS_proportion_LAUR<- ggplot(data=LAUR_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat=\identity\, position = \fill\)+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title=\Laurel Hollow Proportion of body condition scores - 2023\, x =\month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))
  #facet_wrap(~ year)
 #scale_fill_brewer() +
  #facet_wrap(~ site)

BCS_proportion_LAUR
  
pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/LAUR_BCS_proportion.pdf\), height = 7, width = 13)
print(BCS_proportion_LAUR)
dev.off()

<!-- rnb-output-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dVhHNU1RVlZTWDBKRFUxOXdjbTl3YjNKMGFXOXVjend0SUdSbVgweEJWVklnSlQ0bFhHNGdJRzExZEdGMFpTaGpiMjVrYVhScGIyNWZjMk52Y21WZmJuVnRaWEpwWXlBOUlHRnpMbTUxYldWeWFXTW9ZMjl1WkdsMGFXOXVYM05qYjNKbEtTeGpiMjVrYVhScGIyNWZjMk52Y21WZlltbHVJRDBnWTJGelpWOTNhR1Z1S0dOdmJtUnBkR2x2Ymw5elkyOXlaU0E5UFNBeElINGdYRnd4WEZ3c0lHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQXlJSDVjWERKY1hDeGpiMjVrYVhScGIyNWZjMk52Y21VZ1BUMGdNeUIrSUZ4Y00xeGNMQ0JqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTkNCK1hGdzBYRndzWTI5dVpHbDBhVzl1WDNOamIzSmxJRDA5SURVZ2ZpQmNYRFZjWEN3Z1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlEWWdmbHhjTmx4Y0xHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQTNJSDRnWEZ3M1hGd3NJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0E0SUg1Y1hEaGNYQ3dnWTI5dVpHbDBhVzl1WDNOamIzSmxJRDA5T1NCK1hGdzVYRndzSUZSU1ZVVWdmaUJoY3k1amFHRnlZV04wWlhJb1kyOXVaR2wwYVc5dVgzTmpiM0psS1NrcElDVStKVnh1SUNCbmNtOTFjRjlpZVNoemFYUmxMQ0J0YjI1MGFDd2dZMjl1WkdsMGFXOXVYM05qYjNKbFgySnBiaXdnZVdWaGNpa2dKVDRsSUdSd2JIbHlPanB6ZFcxdFlYSnBjMlVvUTI5MWJuUTlJRzRvS1NrZ0pUNGxYRzRnSUhWdVozSnZkWEFvS1NBbFBpVmNiaUFnYlhWMFlYUmxLRkJ5YjNCdmNuUnBiMjRnUFNCRGIzVnVkQzl6ZFcwb1EyOTFiblFwS1Z4dVhHNU1RVlZTWDBKRFUxOXdjbTl3YjNKMGFXOXVjend0SUc1aExtOXRhWFFvVEVGVlVsOUNRMU5mY0hKdmNHOXlkR2x2Ym5NcFhHNWNia3hCVlZKZlFrTlRYM0J5YjNCdmNuUnBiMjV6WEc1Y2JrSkRVMTl3Y205d2IzSjBhVzl1WDB4QlZWSThMU0JuWjNCc2IzUW9aR0YwWVQxTVFWVlNYMEpEVTE5d2NtOXdiM0owYVc5dWN5d2dZV1Z6S0hnOWJXOXVkR2dzSUhrOUlGQnliM0J2Y25ScGIyNHNJR1pwYkd3OVkyOXVaR2wwYVc5dVgzTmpiM0psWDJKcGJpa3BJQ3RjYmlBZ1oyVnZiVjlpWVhJb2QybGtkR2dnUFNBdU5Td2djM1JoZEQxY1hHbGtaVzUwYVhSNVhGd3NJSEJ2YzJsMGFXOXVJRDBnWEZ4bWFXeHNYRndwS3lBZ1hHNGdJQ0IwYUdWdFpWOWlkeWdwSUNzZ0lIUm9aVzFsS0hCaGJtVnNMbWR5YVdRdWJXRnFiM0lnUFNCbGJHVnRaVzUwWDJKc1lXNXJLQ2tzSUhCaGJtVnNMbWR5YVdRdWJXbHViM0lnUFNCbGJHVnRaVzUwWDJKc1lXNXJLQ2twS3lCY2JpQWdkR2hsYldVb1lYaHBjeTUwWlhoMExuZ2dQU0JsYkdWdFpXNTBYM1JsZUhRb1lXNW5iR1VnUFNBME5Td2dkbXAxYzNRZ1BTQXhMQ0JvYW5WemREMHhLU2tyWEc0Z0lHeGhZbk1vZEdsMGJHVTlYRnhNWVhWeVpXd2dTRzlzYkc5M0lGQnliM0J2Y25ScGIyNGdiMllnWW05a2VTQmpiMjVrYVhScGIyNGdjMk52Y21WeklDMGdNakF5TTF4Y0xDQjRJRDFjWEcxdmJuUm9YRndzSUhrZ1BTQmNYRkJ5YjNCdmNuUnBiMjRnYjJZZ1ltOWtlU0JqYjI1a2FYUnBiMjRnYzJOdmNtVnpJRnhjS1NzZ2RHaGxiV1VvWVhocGN5NTBhWFJzWlM1NUlEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0J5Wld3b01TNHpLU3dnWVc1bmJHVWdQVGt3S1N3Z1lYaHBjeTUwYVhSc1pTNTRJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFNBd0tTa3JYRzRnSUhSb1pXMWxLR0Y0YVhNdWRHVjRkRDFsYkdWdFpXNTBYM1JsZUhRb2MybDZaVDB4TWlrcFhHNGdJQ05tWVdObGRGOTNjbUZ3S0g0Z2VXVmhjaWxjYmlBamMyTmhiR1ZmWm1sc2JGOWljbVYzWlhJb0tTQXJYRzRnSUNObVlXTmxkRjkzY21Gd0tINGdjMmwwWlNsY2JseHVRa05UWDNCeWIzQnZjblJwYjI1ZlRFRlZVbHh1SUNCY2JuQmtaaWh3WVhOMFpUQW9jR0YwYUNBOUlGeGNUR0ZpWDBSaGRHRmZWR2x6YzNWbFVISnZZMlZ6YzJsdVp5OXZkWFJ3ZFhSY1hDQXNYRnd2VEVGVlVsOUNRMU5mY0hKdmNHOXlkR2x2Ymk1d1pHWmNYQ2tzSUdobGFXZG9kQ0E5SURjc0lIZHBaSFJvSUQwZ01UTXBYRzV3Y21sdWRDaENRMU5mY0hKdmNHOXlkR2x2Ymw5TVFWVlNLVnh1WkdWMkxtOW1aaWdwWEc1Z1lHQmNibUJnWUNKOSAtLT5cblxuYGBgclxuYGBgclxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIGRmX0xBVVIgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbiwgeWVhcikgJT4lIGRwbHlyOjpzdW1tYXJpc2UoQ291bnQ9IG4oKSkgJT4lXG4gIHVuZ3JvdXAoKSAlPiVcbiAgbXV0YXRlKFByb3BvcnRpb24gPSBDb3VudC9zdW0oQ291bnQpKVxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoTEFVUl9CQ1NfcHJvcG9ydGlvbnMpXG5cbkxBVVJfQkNTX3Byb3BvcnRpb25zXG5cbkJDU19wcm9wb3J0aW9uX0xBVVI8LSBnZ3Bsb3QoZGF0YT1MQVVSX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpKyAgXG4gICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgdmp1c3QgPSAxLCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIFxcKSsgdGhlbWUoYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPTkwKSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSkrXG4gIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMikpXG4gICNmYWNldF93cmFwKH4geWVhcilcbiAjc2NhbGVfZmlsbF9icmV3ZXIoKSArXG4gICNmYWNldF93cmFwKH4gc2l0ZSlcblxuQkNTX3Byb3BvcnRpb25fTEFVUlxuICBcbnBkZihwYXN0ZTAocGF0aCA9IFxcTGFiX0RhdGFfVGlzc3VlUHJvY2Vzc2luZy9vdXRwdXRcXCAsXFwvTEFVUl9CQ1NfcHJvcG9ydGlvbi5wZGZcXCksIGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChCQ1NfcHJvcG9ydGlvbl9MQVVSKVxuZGV2Lm9mZigpXG5gYGBcbmBgYFxuXG48IS0tIHJuYi1zb3VyY2UtZW5kIC0tPlxuIn0= -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIGRmX0xBVVIgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbiwgeWVhcikgJT4lIGRwbHlyOjpzdW1tYXJpc2UoQ291bnQ9IG4oKSkgJT4lXG4gIHVuZ3JvdXAoKSAlPiVcbiAgbXV0YXRlKFByb3BvcnRpb24gPSBDb3VudC9zdW0oQ291bnQpKVxuXG5MQVVSX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoTEFVUl9CQ1NfcHJvcG9ydGlvbnMpXG5cbkxBVVJfQkNTX3Byb3BvcnRpb25zXG5cbkJDU19wcm9wb3J0aW9uX0xBVVI8LSBnZ3Bsb3QoZGF0YT1MQVVSX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpKyAgXG4gICB0aGVtZV9idygpICsgIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCkpKyBcbiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgdmp1c3QgPSAxLCBoanVzdD0xKSkrXG4gIGxhYnModGl0bGU9XFxMYXVyZWwgSG9sbG93IFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIC0gMjAyM1xcLCB4ID1cXG1vbnRoXFwsIHkgPSBcXFByb3BvcnRpb24gb2YgYm9keSBjb25kaXRpb24gc2NvcmVzIFxcKSsgdGhlbWUoYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPTkwKSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSkrXG4gIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMikpXG4gICNmYWNldF93cmFwKH4geWVhcilcbiAjc2NhbGVfZmlsbF9icmV3ZXIoKSArXG4gICNmYWNldF93cmFwKH4gc2l0ZSlcblxuQkNTX3Byb3BvcnRpb25fTEFVUlxuICBcbnBkZihwYXN0ZTAocGF0aCA9IFxcTGFiX0RhdGFfVGlzc3VlUHJvY2Vzc2luZy9vdXRwdXRcXCAsXFwvTEFVUl9CQ1NfcHJvcG9ydGlvbi5wZGZcXCksIGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChCQ1NfcHJvcG9ydGlvbl9MQVVSKVxuZGV2Lm9mZigpXG5gYGBcbmBgYCJ9 -->

```r
```r

LAUR_BCS_proportions<- df_LAUR %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

LAUR_BCS_proportions<- na.omit(LAUR_BCS_proportions)

LAUR_BCS_proportions

BCS_proportion_LAUR<- ggplot(data=LAUR_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat=\identity\, position = \fill\)+  
   theme_bw() +  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
  labs(title=\Laurel Hollow Proportion of body condition scores - 2023\, x =\month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
  theme(axis.text=element_text(size=12))
  #facet_wrap(~ year)
 #scale_fill_brewer() +
  #facet_wrap(~ site)

BCS_proportion_LAUR
  
pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/LAUR_BCS_proportion.pdf\), height = 7, width = 13)
print(BCS_proportion_LAUR)
dev.off()
```
```

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-text-begin -->

# % of scores >3 at LAUR june - november 

<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVYRzRqY21WdGIzWnBibWNnTURneU0wWkZUa05mTWpnZ1ltTnpJSGRoY3lCT1FWeHVaR1pmVEVGVlVpNGdQQzBnWkdaZlRFRlZVaVUrSldacGJIUmxjaWdoY205M1gyNTFiV0psY2lncElDVnBiaVVnWXlneUtTbGNibHh1VEVGVlVsOXdjbVZqWlc1MFgyZHlaV0YwWlhKZk15QThMU0JrWmw5TVFWVlNMaUFsUGlWY2JpQWdaSEJzZVhJNk9tZHliM1Z3WDJKNUtHMXZiblJvTENCemFYUmxMQ0I1WldGeUtTQWxQaVZjYmlBZ1pIQnNlWEk2T25OMWJXMWhjbWx6WlNoUVpYSmpaVzUwWVdkbElEMGdiV1ZoYmloamIyNWthWFJwYjI1ZmMyTnZjbVVnUEQwZ015a3FNVEF3S1Z4dVRFRlZVbDl3Y21WalpXNTBYMmR5WldGMFpYSmZNMXh1WEc1TVFWVlNYMEpEVTE5d1pYSmpaVzUwWVdkbElEd3RJRXhCVlZKZmNISmxZMlZ1ZEY5bmNtVmhkR1Z5WHpNZ0pUNGxYRzRnSUdkbmNHeHZkQ2hoWlhNb2VDQTlJRzF2Ym5Sb0xDQjVJRDBnVUdWeVkyVnVkR0ZuWlNrcElDdGNiaUFnSTJkbGIyMWZZbUZ5S0hkcFpIUm9JRDBnTUM0MUxDQnpkR0YwSUQwZ1hDSnBaR1Z1ZEdsMGVWd2lMQ0J3YjNOcGRHbHZiaUE5SUZ3aVptbHNiRndpS1NBclhHNGdJQ0FnWjJWdmJWOWpiMndvS1N0Y2JpQWdkR2hsYldWZlluY29LU0FyWEc0Z0lIUm9aVzFsS0Z4dUlDQWdJSEJoYm1Wc0xtZHlhV1F1YldGcWIzSWdQU0JsYkdWdFpXNTBYMkpzWVc1cktDa3NYRzRnSUNBZ2NHRnVaV3d1WjNKcFpDNXRhVzV2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1N4Y2JpQWdJQ0JoZUdsekxuUmxlSFF1ZUNBOUlHVnNaVzFsYm5SZmRHVjRkQ2hoYm1kc1pTQTlJRFExTENCMmFuVnpkQ0E5SURFc0lHaHFkWE4wSUQwZ01Ta3NYRzRnSUNBZ1lYaHBjeTUwYVhSc1pTNTVJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFNBNU1Da3NYRzRnSUNBZ1lYaHBjeTUwYVhSc1pTNTRJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNCeVpXd29NUzR6S1N3Z1lXNW5iR1VnUFNBd0tTeGNiaUFnSUNCaGVHbHpMblJsZUhRZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvYzJsNlpTQTlJREV5S1Z4dUlDQXBJQ3RjYmlBZ2JHRmljeWhjYmlBZ0lDQjBhWFJzWlNBOUlGd2lSMjlzWkNCVGRHRnlJRUpsWVdOb0lDVWdiMllnWW05a2VTQmpiMjVrYVhScGIyNGdjMk52Y21WeklENDlJRE5jSWl4Y2JpQWdJQ0I0SUQwZ1hDSnRiMjUwYUZ3aUxGeHVJQ0FnSUhrZ1BTQmNJbEJsY21ObGJuUmhaMlVnYjJZZ1ltOWtlU0JqYjI1a2FYUnBiMjRnYzJOdmNtVnpJRDQ5SUROY0lseHVJQ0FwSUNNclptRmpaWFJmZDNKaGNDaCtlV1ZoY2lsY2JpQWdJeXNnYzJOaGJHVmZlRjlqYjI1MGFXNTFiM1Z6S0NCaWNtVmhhM01nUFNCelpYRW9OU3d4TWl3Z1lua2dQVEVwSUNrclhHNGdJQ041YkdsdEtEQXNOakFwS3lCY2JpQWdJM05qWVd4bFgzbGZZMjl1ZEdsdWRXOTFjeWhzYVcxcGRITWdQU0JqS0RBc01UQXdLU3dnWW5KbFlXdHpJRDBnYzJWeEtEQXNNVEF3TENCaWVTQTlJREV3S1NrZ1hHNWNiaUFnSUNCY2JreEJWVkpmUWtOVFgzQmxjbU5sYm5SaFoyVmNibHh1Y0dSbUtIQmhjM1JsTUNod1lYUm9JRDBnWENKTVlXSmZSR0YwWVY5VWFYTnpkV1ZRY205alpYTnphVzVuTDI5MWRIQjFkRndpTENCY0lpOU1RVlZTWDBKRFUxOXdaWEpqWlc1MFlXZGxMbkJrWmx3aUtTeG9aV2xuYUhRZ1BTQTNMQ0IzYVdSMGFDQTlJREV6S1Z4dWNISnBiblFvVEVGVlVsOUNRMU5mY0dWeVkyVnVkR0ZuWlNsY2JtUmxkaTV2Wm1Zb0tTQmNibHh1WEc1Y2JtQmdZQ0o5IC0tPlxuXG5gYGByXG5cbiNyZW1vdmluZyAwODIzRkVOQ18yOCBiY3Mgd2FzIE5BXG5kZl9MQVVSLiA8LSBkZl9MQVVSJT4lZmlsdGVyKCFyb3dfbnVtYmVyKCkgJWluJSBjKDIpKVxuXG5MQVVSX3ByZWNlbnRfZ3JlYXRlcl8zIDwtIGRmX0xBVVIuICU+JVxuICBkcGx5cjo6Z3JvdXBfYnkobW9udGgsIHNpdGUsIHllYXIpICU+JVxuICBkcGx5cjo6c3VtbWFyaXNlKFBlcmNlbnRhZ2UgPSBtZWFuKGNvbmRpdGlvbl9zY29yZSA8PSAzKSoxMDApXG5MQVVSX3ByZWNlbnRfZ3JlYXRlcl8zXG5cbkxBVVJfQkNTX3BlcmNlbnRhZ2UgPC0gTEFVUl9wcmVjZW50X2dyZWF0ZXJfMyAlPiVcbiAgZ2dwbG90KGFlcyh4ID0gbW9udGgsIHkgPSBQZXJjZW50YWdlKSkgK1xuICAjZ2VvbV9iYXIod2lkdGggPSAwLjUsIHN0YXQgPSBcXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwpICtcbiAgICBnZW9tX2NvbCgpK1xuICB0aGVtZV9idygpICtcbiAgdGhlbWUoXG4gICAgcGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSxcbiAgICBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpLFxuICAgIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIHZqdXN0ID0gMSwgaGp1c3QgPSAxKSxcbiAgICBheGlzLnRpdGxlLnkgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjMpLCBhbmdsZSA9IDkwKSxcbiAgICBheGlzLnRpdGxlLnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjMpLCBhbmdsZSA9IDApLFxuICAgIGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTIpXG4gICkgK1xuICBsYWJzKFxuICAgIHRpdGxlID0gXFxHb2xkIFN0YXIgQmVhY2ggJSBvZiBib2R5IGNvbmRpdGlvbiBzY29yZXMgPj0gM1xcLFxuICAgIHggPSBcXG1vbnRoXFwsXG4gICAgeSA9IFxcUGVyY2VudGFnZSBvZiBib2R5IGNvbmRpdGlvbiBzY29yZXMgPj0gM1xcXG4gICkgIytmYWNldF93cmFwKH55ZWFyKVxuICAjKyBzY2FsZV94X2NvbnRpbnVvdXMoIGJyZWFrcyA9IHNlcSg1LDEyLCBieSA9MSkgKStcbiAgI3lsaW0oMCw2MCkrIFxuICAjc2NhbGVfeV9jb250aW51b3VzKGxpbWl0cyA9IGMoMCwxMDApLCBicmVha3MgPSBzZXEoMCwxMDAsIGJ5ID0gMTApKSBcblxuICAgIFxuTEFVUl9CQ1NfcGVyY2VudGFnZVxuXG5wZGYocGFzdGUwKHBhdGggPSBcXExhYl9EYXRhX1Rpc3N1ZVByb2Nlc3Npbmcvb3V0cHV0XFwsIFxcL0xBVVJfQkNTX3BlcmNlbnRhZ2UucGRmXFwpLGhlaWdodCA9IDcsIHdpZHRoID0gMTMpXG5wcmludChMQVVSX0JDU19wZXJjZW50YWdlKVxuZGV2Lm9mZigpIFxuXG5cblxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->

#removing 0823FENC_28 bcs was NA
df_LAUR. <- df_LAUR%>%filter(!row_number() %in% c(2))

LAUR_precent_greater_3 <- df_LAUR. %>%
  dplyr::group_by(month, site, year) %>%
  dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
LAUR_precent_greater_3

LAUR_BCS_percentage <- LAUR_precent_greater_3 %>%
  ggplot(aes(x = month, y = Percentage)) +
  #geom_bar(width = 0.5, stat = \identity\, position = \fill\) +
    geom_col()+
  theme_bw() +
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
    axis.title.y = element_text(size = rel(1.3), angle = 90),
    axis.title.x = element_text(size = rel(1.3), angle = 0),
    axis.text = element_text(size = 12)
  ) +
  labs(
    title = \Gold Star Beach % of body condition scores >= 3\,
    x = \month\,
    y = \Percentage of body condition scores >= 3\
  ) #+facet_wrap(~year)
  #+ scale_x_continuous( breaks = seq(5,12, by =1) )+
  #ylim(0,60)+ 
  #scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10)) 

    
LAUR_BCS_percentage

pdf(paste0(path = \Lab_Data_TissueProcessing/output\, \/LAUR_BCS_percentage.pdf\),height = 7, width = 13)
print(LAUR_BCS_percentage)
dev.off() 

<!-- rnb-output-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dVhHNGpjbVZ0YjNacGJtY2dNRGd5TTBaRlRrTmZNamdnWW1OeklIZGhjeUJPUVZ4dVpHWmZURUZWVWk0Z1BDMGdaR1pmVEVGVlVpVStKV1pwYkhSbGNpZ2hjbTkzWDI1MWJXSmxjaWdwSUNWcGJpVWdZeWd5S1NsY2JseHVURUZWVWw5d2NtVmpaVzUwWDJkeVpXRjBaWEpmTXlBOExTQmtabDlNUVZWU0xpQWxQaVZjYmlBZ1pIQnNlWEk2T21keWIzVndYMko1S0cxdmJuUm9MQ0J6YVhSbExDQjVaV0Z5S1NBbFBpVmNiaUFnWkhCc2VYSTZPbk4xYlcxaGNtbHpaU2hRWlhKalpXNTBZV2RsSUQwZ2JXVmhiaWhqYjI1a2FYUnBiMjVmYzJOdmNtVWdQRDBnTXlrcU1UQXdLVnh1VEVGVlVsOXdjbVZqWlc1MFgyZHlaV0YwWlhKZk0xeHVYRzVNUVZWU1gwSkRVMTl3WlhKalpXNTBZV2RsSUR3dElFeEJWVkpmY0hKbFkyVnVkRjluY21WaGRHVnlYek1nSlQ0bFhHNGdJR2RuY0d4dmRDaGhaWE1vZUNBOUlHMXZiblJvTENCNUlEMGdVR1Z5WTJWdWRHRm5aU2twSUN0Y2JpQWdJMmRsYjIxZlltRnlLSGRwWkhSb0lEMGdNQzQxTENCemRHRjBJRDBnWEZ4cFpHVnVkR2wwZVZ4Y0xDQndiM05wZEdsdmJpQTlJRnhjWm1sc2JGeGNLU0FyWEc0Z0lDQWdaMlZ2YlY5amIyd29LU3RjYmlBZ2RHaGxiV1ZmWW5jb0tTQXJYRzRnSUhSb1pXMWxLRnh1SUNBZ0lIQmhibVZzTG1keWFXUXViV0ZxYjNJZ1BTQmxiR1Z0Wlc1MFgySnNZVzVyS0Nrc1hHNGdJQ0FnY0dGdVpXd3VaM0pwWkM1dGFXNXZjaUE5SUdWc1pXMWxiblJmWW14aGJtc29LU3hjYmlBZ0lDQmhlR2x6TG5SbGVIUXVlQ0E5SUdWc1pXMWxiblJmZEdWNGRDaGhibWRzWlNBOUlEUTFMQ0IyYW5WemRDQTlJREVzSUdocWRYTjBJRDBnTVNrc1hHNGdJQ0FnWVhocGN5NTBhWFJzWlM1NUlEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0J5Wld3b01TNHpLU3dnWVc1bmJHVWdQU0E1TUNrc1hHNGdJQ0FnWVhocGN5NTBhWFJzWlM1NElEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0J5Wld3b01TNHpLU3dnWVc1bmJHVWdQU0F3S1N4Y2JpQWdJQ0JoZUdsekxuUmxlSFFnUFNCbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlNBOUlERXlLVnh1SUNBcElDdGNiaUFnYkdGaWN5aGNiaUFnSUNCMGFYUnNaU0E5SUZ4Y1IyOXNaQ0JUZEdGeUlFSmxZV05vSUNVZ2IyWWdZbTlrZVNCamIyNWthWFJwYjI0Z2MyTnZjbVZ6SUQ0OUlETmNYQ3hjYmlBZ0lDQjRJRDBnWEZ4dGIyNTBhRnhjTEZ4dUlDQWdJSGtnUFNCY1hGQmxjbU5sYm5SaFoyVWdiMllnWW05a2VTQmpiMjVrYVhScGIyNGdjMk52Y21WeklENDlJRE5jWEZ4dUlDQXBJQ01yWm1GalpYUmZkM0poY0NoK2VXVmhjaWxjYmlBZ0l5c2djMk5oYkdWZmVGOWpiMjUwYVc1MWIzVnpLQ0JpY21WaGEzTWdQU0J6WlhFb05Td3hNaXdnWW5rZ1BURXBJQ2tyWEc0Z0lDTjViR2x0S0RBc05qQXBLeUJjYmlBZ0kzTmpZV3hsWDNsZlkyOXVkR2x1ZFc5MWN5aHNhVzFwZEhNZ1BTQmpLREFzTVRBd0tTd2dZbkpsWVd0eklEMGdjMlZ4S0RBc01UQXdMQ0JpZVNBOUlERXdLU2tnWEc1Y2JpQWdJQ0JjYmt4QlZWSmZRa05UWDNCbGNtTmxiblJoWjJWY2JseHVjR1JtS0hCaGMzUmxNQ2h3WVhSb0lEMGdYRnhNWVdKZlJHRjBZVjlVYVhOemRXVlFjbTlqWlhOemFXNW5MMjkxZEhCMWRGeGNMQ0JjWEM5TVFWVlNYMEpEVTE5d1pYSmpaVzUwWVdkbExuQmtabHhjS1N4b1pXbG5hSFFnUFNBM0xDQjNhV1IwYUNBOUlERXpLVnh1Y0hKcGJuUW9URUZWVWw5Q1ExTmZjR1Z5WTJWdWRHRm5aU2xjYm1SbGRpNXZabVlvS1NCY2JseHVYRzVjYm1CZ1lGeHVZR0JnSW4wPSAtLT5cblxuYGBgclxuYGBgclxuXG4jcmVtb3ZpbmcgMDgyM0ZFTkNfMjggYmNzIHdhcyBOQVxuZGZfTEFVUi4gPC0gZGZfTEFVUiU+JWZpbHRlcighcm93X251bWJlcigpICVpbiUgYygyKSlcblxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfMyA8LSBkZl9MQVVSLiAlPiVcbiAgZHBseXI6Omdyb3VwX2J5KG1vbnRoLCBzaXRlLCB5ZWFyKSAlPiVcbiAgZHBseXI6OnN1bW1hcmlzZShQZXJjZW50YWdlID0gbWVhbihjb25kaXRpb25fc2NvcmUgPD0gMykqMTAwKVxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfM1xuXG5MQVVSX0JDU19wZXJjZW50YWdlIDwtIExBVVJfcHJlY2VudF9ncmVhdGVyXzMgJT4lXG4gIGdncGxvdChhZXMoeCA9IG1vbnRoLCB5ID0gUGVyY2VudGFnZSkpICtcbiAgI2dlb21fYmFyKHdpZHRoID0gMC41LCBzdGF0ID0gXFxpZGVudGl0eVxcLCBwb3NpdGlvbiA9IFxcZmlsbFxcKSArXG4gICAgZ2VvbV9jb2woKStcbiAgdGhlbWVfYncoKSArXG4gIHRoZW1lKFxuICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksXG4gICAgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSxcbiAgICBheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCB2anVzdCA9IDEsIGhqdXN0ID0gMSksXG4gICAgYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSA5MCksXG4gICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSxcbiAgICBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDEyKVxuICApICtcbiAgbGFicyhcbiAgICB0aXRsZSA9IFxcR29sZCBTdGFyIEJlYWNoICUgb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXCxcbiAgICB4ID0gXFxtb250aFxcLFxuICAgIHkgPSBcXFBlcmNlbnRhZ2Ugb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXFxuICApICMrZmFjZXRfd3JhcCh+eWVhcilcbiAgIysgc2NhbGVfeF9jb250aW51b3VzKCBicmVha3MgPSBzZXEoNSwxMiwgYnkgPTEpICkrXG4gICN5bGltKDAsNjApKyBcbiAgI3NjYWxlX3lfY29udGludW91cyhsaW1pdHMgPSBjKDAsMTAwKSwgYnJlYWtzID0gc2VxKDAsMTAwLCBieSA9IDEwKSkgXG5cbiAgICBcbkxBVVJfQkNTX3BlcmNlbnRhZ2VcblxucGRmKHBhc3RlMChwYXRoID0gXFxMYWJfRGF0YV9UaXNzdWVQcm9jZXNzaW5nL291dHB1dFxcLCBcXC9MQVVSX0JDU19wZXJjZW50YWdlLnBkZlxcKSxoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxucHJpbnQoTEFVUl9CQ1NfcGVyY2VudGFnZSlcbmRldi5vZmYoKSBcblxuXG5cbmBgYFxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuXG4jcmVtb3ZpbmcgMDgyM0ZFTkNfMjggYmNzIHdhcyBOQVxuZGZfTEFVUi4gPC0gZGZfTEFVUiU+JWZpbHRlcighcm93X251bWJlcigpICVpbiUgYygyKSlcblxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfMyA8LSBkZl9MQVVSLiAlPiVcbiAgZHBseXI6Omdyb3VwX2J5KG1vbnRoLCBzaXRlLCB5ZWFyKSAlPiVcbiAgZHBseXI6OnN1bW1hcmlzZShQZXJjZW50YWdlID0gbWVhbihjb25kaXRpb25fc2NvcmUgPD0gMykqMTAwKVxuTEFVUl9wcmVjZW50X2dyZWF0ZXJfM1xuXG5MQVVSX0JDU19wZXJjZW50YWdlIDwtIExBVVJfcHJlY2VudF9ncmVhdGVyXzMgJT4lXG4gIGdncGxvdChhZXMoeCA9IG1vbnRoLCB5ID0gUGVyY2VudGFnZSkpICtcbiAgI2dlb21fYmFyKHdpZHRoID0gMC41LCBzdGF0ID0gXFxpZGVudGl0eVxcLCBwb3NpdGlvbiA9IFxcZmlsbFxcKSArXG4gICAgZ2VvbV9jb2woKStcbiAgdGhlbWVfYncoKSArXG4gIHRoZW1lKFxuICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksXG4gICAgcGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSxcbiAgICBheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCB2anVzdCA9IDEsIGhqdXN0ID0gMSksXG4gICAgYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSA5MCksXG4gICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSByZWwoMS4zKSwgYW5nbGUgPSAwKSxcbiAgICBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDEyKVxuICApICtcbiAgbGFicyhcbiAgICB0aXRsZSA9IFxcR29sZCBTdGFyIEJlYWNoICUgb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXCxcbiAgICB4ID0gXFxtb250aFxcLFxuICAgIHkgPSBcXFBlcmNlbnRhZ2Ugb2YgYm9keSBjb25kaXRpb24gc2NvcmVzID49IDNcXFxuICApICMrZmFjZXRfd3JhcCh+eWVhcilcbiAgIysgc2NhbGVfeF9jb250aW51b3VzKCBicmVha3MgPSBzZXEoNSwxMiwgYnkgPTEpICkrXG4gICN5bGltKDAsNjApKyBcbiAgI3NjYWxlX3lfY29udGludW91cyhsaW1pdHMgPSBjKDAsMTAwKSwgYnJlYWtzID0gc2VxKDAsMTAwLCBieSA9IDEwKSkgXG5cbiAgICBcbkxBVVJfQkNTX3BlcmNlbnRhZ2VcblxucGRmKHBhc3RlMChwYXRoID0gXFxMYWJfRGF0YV9UaXNzdWVQcm9jZXNzaW5nL291dHB1dFxcLCBcXC9MQVVSX0JDU19wZXJjZW50YWdlLnBkZlxcKSxoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxucHJpbnQoTEFVUl9CQ1NfcGVyY2VudGFnZSlcbmRldi5vZmYoKSBcblxuXG5cbmBgYFxuYGBgIn0= -->

```r
```r

#removing 0823FENC_28 bcs was NA
df_LAUR. <- df_LAUR%>%filter(!row_number() %in% c(2))

LAUR_precent_greater_3 <- df_LAUR. %>%
  dplyr::group_by(month, site, year) %>%
  dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
LAUR_precent_greater_3

LAUR_BCS_percentage <- LAUR_precent_greater_3 %>%
  ggplot(aes(x = month, y = Percentage)) +
  #geom_bar(width = 0.5, stat = \identity\, position = \fill\) +
    geom_col()+
  theme_bw() +
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
    axis.title.y = element_text(size = rel(1.3), angle = 90),
    axis.title.x = element_text(size = rel(1.3), angle = 0),
    axis.text = element_text(size = 12)
  ) +
  labs(
    title = \Gold Star Beach % of body condition scores >= 3\,
    x = \month\,
    y = \Percentage of body condition scores >= 3\
  ) #+facet_wrap(~year)
  #+ scale_x_continuous( breaks = seq(5,12, by =1) )+
  #ylim(0,60)+ 
  #scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10)) 

    
LAUR_BCS_percentage

pdf(paste0(path = \Lab_Data_TissueProcessing/output\, \/LAUR_BCS_percentage.pdf\),height = 7, width = 13)
print(LAUR_BCS_percentage)
dev.off() 

```
```

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-text-begin -->



<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVJME52WkdsdVp5Qm5jbUZ3YUhNZ1ptOXlJRXhKVXlCRGIyNW1aWEpsYm1ObElDMGdiMjVzZVNBeU1ESXpJR1JoZEdGY2JtUmhkR0ZmTWpBeU15QThMU0JrWVhSaFgyRnNiQ1UrSlNBZ1ptbHNkR1Z5S0NGNVpXRnlQVDFjSWpJd01qUmNJaWxjYmx4dVpHWmZRa05UWDNCeWIzQnZjblJwYjI1elBDMGdaR0YwWVY4eU1ESXpJQ1UrSlZ4dUlDQnRkWFJoZEdVb1kyOXVaR2wwYVc5dVgzTmpiM0psWDI1MWJXVnlhV01nUFNCaGN5NXVkVzFsY21saktHTnZibVJwZEdsdmJsOXpZMjl5WlNrc1kyOXVaR2wwYVc5dVgzTmpiM0psWDJKcGJpQTlJR05oYzJWZmQyaGxiaWhqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTVNCK0lGd2lNVndpTENCamIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ01pQitYQ0l5WENJc1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlETWdmaUJjSWpOY0lpd2dZMjl1WkdsMGFXOXVYM05qYjNKbElEMDlJRFFnZmx3aU5Gd2lMR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0ExSUg0Z1hDSTFYQ0lzSUdOdmJtUnBkR2x2Ymw5elkyOXlaU0E5UFNBMklINWNJalpjSWl4amIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ055QitJRndpTjF3aUxDQmpiMjVrYVhScGIyNWZjMk52Y21VZ1BUMGdPQ0IrWENJNFhDSXNJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQVGtnZmx3aU9Wd2lMQ0JVVWxWRklINGdZWE11WTJoaGNtRmpkR1Z5S0dOdmJtUnBkR2x2Ymw5elkyOXlaU2twS1NBbFBpVmNiaUFnWjNKdmRYQmZZbmtvYzJsMFpTd2diVzl1ZEdnc0lHTnZibVJwZEdsdmJsOXpZMjl5WlY5aWFXNHBJQ1UrSlNCa2NHeDVjaUE2T25OMWJXMWhjbWx6WlNoRGIzVnVkRDBnYmlncEtTQWxQaVZjYmlBZ2RXNW5jbTkxY0NncElDVStKVnh1SUNCdGRYUmhkR1VvVUhKdmNHOXlkR2x2YmlBOUlFTnZkVzUwTDNOMWJTaERiM1Z1ZENrcFhHNWNibVJtWDBKRFUxOXdjbTl3YjNKMGFXOXVjend0SUc1aExtOXRhWFFvWkdaZlFrTlRYM0J5YjNCdmNuUnBiMjV6S1Z4dVhHNWtabDlDUTFOZmNISnZjRzl5ZEdsdmJuTWtiVzl1ZEdnZ1BDMGdabUZqZEc5eUtHUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN5UnRiMjUwYUN3Z2JHVjJaV3h6SUQwZ1l5aGNJalZjSWl4Y0lqWmNJaXhjSWpkY0lpd2dYQ0k0WENJc0lGd2lPVndpTENCY0lqRXdYQ0lzSUZ3aU1URmNJaWtzWEc0Z0lDQWdJQ0FnSUd4aFltVnNjejFqS0Z3aVRXRjVYQ0lzSUZ3aVNuVnVaVndpTENCY0lrcDFiSGxjSWl3Z1hDSkJkV2RjSWl3Z1hDSlRaWEIwWENJc0lGd2lUMk4wWENJc0lGd2lUbTkyWENJcEtWeHVaR1pmUWtOVFgzQnliM0J2Y25ScGIyNXpKSE5wZEdVZ1BDMGdabUZqZEc5eUtHUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN5UnphWFJsTENCc1pYWmxiSE1nUFNCaktGd2lRVk5JUTF3aUxGd2lSa1ZPUTF3aUxGd2lSMDlNUkZ3aUxDQmNJa3hCVlZKY0lpa3NYRzRnSUNBZ0lDQWdJR3hoWW1Wc2N6MWpLRndpUVhOb0lFTnlaV1ZyWENJc0lGd2lSbVZ1WTJVZ1EzSmxaV3RjSWl3Z1hDSkhiMnhrSUZOMFlYSmNJaXdnWENKTVlYVnlaV3dnU0c5c2JHOTNYQ0lwS1Z4dVhHNWlZM011Y0hKdmNDNHlNREl6SUR3dElHZG5jR3h2ZENoa1lYUmhQV1JtWDBKRFUxOXdjbTl3YjNKMGFXOXVjeXdnWVdWektIZzliVzl1ZEdnc0lIazlJRkJ5YjNCdmNuUnBiMjRzSUdacGJHdzlZMjl1WkdsMGFXOXVYM05qYjNKbFgySnBiaWtwSUN0Y2JpQWdaMlZ2YlY5aVlYSW9kMmxrZEdnZ1BTQXVOU3dnYzNSaGREMWNJbWxrWlc1MGFYUjVYQ0lzSUhCdmMybDBhVzl1SUQwZ1hDSm1hV3hzWENJc0lHTnZiRzkxY2lBOUlGd2lZbXhoWTJ0Y0lpa3JJQ0JjYmlBZ0lIUm9aVzFsWDJKM0tDa2dLeUFnZEdobGJXVW9jR0Z1Wld3dVozSnBaQzV0YVc1dmNpQTlJR1ZzWlcxbGJuUmZZbXhoYm1zb0tTa3JJRnh1SUNCMGFHVnRaU2hoZUdsekxuUmxlSFF1ZUNBOUlHVnNaVzFsYm5SZmRHVjRkQ2hoYm1kc1pTQTlJRFExTENCMmFuVnpkQ0E5SURFc0lHaHFkWE4wUFRFc0lITnBlbVVnUFNCeVpXd29MamtwS1NrclhHNGdJSFJvWlcxbEtIQmhibVZzTG1keWFXUXViV0ZxYjNJdWVDQTlaV3hsYldWdWRGOWliR0Z1YXlncEtTdGNiaUFnYkdGaWN5aDBhWFJzWlQxY0lrSnZaSGtnWTI5dVpHbDBhVzl1SUhOamIzSmxjeUJwYmlBeU1ESXpYQ0lzSUhnZ1BWd2lUVzl1ZEdoY0lpd2dlU0E5SUZ3aVVISnZjRzl5ZEdsdmJpQnZaaUJpYjJSNUlHTnZibVJwZEdsdmJpQnpZMjl5WlhNZ1hDSXBLeUIwYUdWdFpTaGhlR2x6TG5ScGRHeGxMbmtnUFNCbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlNBOUlISmxiQ2d4TGpJcExDQmhibWRzWlNBOU9UQXBMQ0JoZUdsekxuUnBkR3hsTG5nZ1BTQmxiR1Z0Wlc1MFgzUmxlSFFvYzJsNlpTQTlJSEpsYkNneExqSXBMQ0JoYm1kc1pTQTlJREFwS1NzZ2MyTmhiR1ZmWm1sc2JGOWljbVYzWlhJb2NHRnNaWFIwWlNBOUlGd2lRbXgxWlhOY0lpa2dLMXh1SUNCMGFHVnRaU2h6ZEhKcGNDNTBaWGgwSUQwZ1pXeGxiV1Z1ZEY5MFpYaDBLSE5wZW1VZ1BTQXhOU2twSzF4dUlDQjBhR1Z0WlNBb2JHVm5aVzVrTG5ScGRHeGxJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNBeE9Da3BLMXh1SUNCMGFHVnRaU2gwYVhSc1pTQTlJR1ZzWlcxbGJuUmZkR1Y0ZENoemFYcGxJRDBnTVRjcEtTdGNiaUFnZEdobGJXVW9ZWGhwY3k1MFpYaDBQV1ZzWlcxbGJuUmZkR1Y0ZENoemFYcGxQVEUzS1NrclhHNGdJR3hoWW5Nb1ptbHNiQ0E5WENKQ2IyUjVJRU52Ym1ScGRHbHZiaUJUWTI5eVpWd2lLU3RjYmlBZ1ptRmpaWFJmZDNKaGNDaCtJSE5wZEdVZ0tWeHVZbU56TG5CeWIzQXVNakF5TTF4dVhHNGpjR1JtS0hCaGMzUmxNQ2h3WVhSb0lEMGdYQ0pNWVdKZlJHRjBZVjlVYVhOemRXVlFjbTlqWlhOemFXNW5MMjkxZEhCMWRGd2lJQ3hjSWk5Q1ExTmZjSEp2Y0c5eWRHbHZibDh5TURJekxuQmtabHdpS1N3Z2FHVnBaMmgwSUQwZ055d2dkMmxrZEdnZ1BTQXhNeWxjYmlOd2NtbHVkQ2hpWTNNdWNISnZjQzR5TURJektWeHVJMlJsZGk1dlptWW9LVnh1WEc0alFrTlRYM0J5YjNCdmNuUnBiMjVmWVd4c0szTmpZV3hsWDJacGJHeGZiV0Z1ZFdGc0tIWmhiSFZsY3lBOUlHTW9YQ0lqTURBelF6TXdYQ0lzWENJak1ERTJOalZGWENJc1hDSWpPREJEUkVNeFhDSXNYQ0lqUXpkRlFVVTFYQ0lzWENJalJqWkZPRU16WENJc1hDSWpSRVpETWpkRVhDSXNYQ0lqUWtZNE1USkVYQ0lzSUZ3aUl6aEROVEV3UVZ3aUxGd2lJelUwTXpBd05Wd2lJQ2twWEc0amMyTmhiR1ZmWm1sc2JGOXRZVzUxWVd3b2RtRnNkV1Z6SUQwZ1l5aGNJaU13TURORE16QmNJaXhjSWlNd01UWTJOVVZjSWl4Y0lpTTRNRU5FUXpGY0lpeGNJaU5ETjBWQlJUVmNJaXhjSW1GNmRYSmxNVndpTEZ3aWMyeGhkR1ZuY21GNU1Wd2lMRndpYzJ4aGRHVm5jbUY1TTF3aUxDQmNJbk5zWVhSbFozSmhlVFJjSWl3Z1hDSm5jbUY1TWpWY0lpQXBLU0JjYmx4dVppQThMU0JtZFc1amRHbHZiaWh3WVd3cElHSnlaWGRsY2k1d1lXd29ZbkpsZDJWeUxuQmhiQzVwYm1adlczQmhiQ3dnWENKdFlYaGpiMnh2Y25OY0lsMHNJSEJoYkNsY2JpaGpiMnh6SUR3dElHWW9YQ0paYkVkdVFuVmNJaWtwWEc1Z1lHQWlmUT09IC0tPlxuXG5gYGByXG4jQ29kaW5nIGdyYXBocyBmb3IgTElTIENvbmZlcmVuY2UgLSBvbmx5IDIwMjMgZGF0YVxuZGF0YV8yMDIzIDwtIGRhdGFfYWxsJT4lICBmaWx0ZXIoIXllYXI9PVxcMjAyNFxcKVxuXG5kZl9CQ1NfcHJvcG9ydGlvbnM8LSBkYXRhXzIwMjMgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbikgJT4lIGRwbHlyIDo6c3VtbWFyaXNlKENvdW50PSBuKCkpICU+JVxuICB1bmdyb3VwKCkgJT4lXG4gIG11dGF0ZShQcm9wb3J0aW9uID0gQ291bnQvc3VtKENvdW50KSlcblxuZGZfQkNTX3Byb3BvcnRpb25zPC0gbmEub21pdChkZl9CQ1NfcHJvcG9ydGlvbnMpXG5cbmRmX0JDU19wcm9wb3J0aW9ucyRtb250aCA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJG1vbnRoLCBsZXZlbHMgPSBjKFxcNVxcLFxcNlxcLFxcN1xcLCBcXDhcXCwgXFw5XFwsIFxcMTBcXCwgXFwxMVxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxNYXlcXCwgXFxKdW5lXFwsIFxcSnVseVxcLCBcXEF1Z1xcLCBcXFNlcHRcXCwgXFxPY3RcXCwgXFxOb3ZcXCkpXG5kZl9CQ1NfcHJvcG9ydGlvbnMkc2l0ZSA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJHNpdGUsIGxldmVscyA9IGMoXFxBU0hDXFwsXFxGRU5DXFwsXFxHT0xEXFwsIFxcTEFVUlxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxBc2ggQ3JlZWtcXCwgXFxGZW5jZSBDcmVla1xcLCBcXEdvbGQgU3RhclxcLCBcXExhdXJlbCBIb2xsb3dcXCkpXG5cbmJjcy5wcm9wLjIwMjMgPC0gZ2dwbG90KGRhdGE9ZGZfQkNTX3Byb3BvcnRpb25zLCBhZXMoeD1tb250aCwgeT0gUHJvcG9ydGlvbiwgZmlsbD1jb25kaXRpb25fc2NvcmVfYmluKSkgK1xuICBnZW9tX2Jhcih3aWR0aCA9IC41LCBzdGF0PVxcaWRlbnRpdHlcXCwgcG9zaXRpb24gPSBcXGZpbGxcXCwgY29sb3VyID0gXFxibGFja1xcKSsgIFxuICAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIHZqdXN0ID0gMSwgaGp1c3Q9MSwgc2l6ZSA9IHJlbCguOSkpKStcbiAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvci54ID1lbGVtZW50X2JsYW5rKCkpK1xuICBsYWJzKHRpdGxlPVxcQm9keSBjb25kaXRpb24gc2NvcmVzIGluIDIwMjNcXCwgeCA9XFxNb250aFxcLCB5ID0gXFxQcm9wb3J0aW9uIG9mIGJvZHkgY29uZGl0aW9uIHNjb3JlcyBcXCkrIHRoZW1lKGF4aXMudGl0bGUueSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID05MCksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID0gMCkpKyBzY2FsZV9maWxsX2JyZXdlcihwYWxldHRlID0gXFxCbHVlc1xcKSArXG4gIHRoZW1lKHN0cmlwLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkrXG4gIHRoZW1lIChsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE4KSkrXG4gIHRoZW1lKHRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNykpK1xuICB0aGVtZShheGlzLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9MTcpKStcbiAgbGFicyhmaWxsID1cXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwpK1xuICBmYWNldF93cmFwKH4gc2l0ZSApXG5iY3MucHJvcC4yMDIzXG5cbiNwZGYocGFzdGUwKHBhdGggPSBcXExhYl9EYXRhX1Rpc3N1ZVByb2Nlc3Npbmcvb3V0cHV0XFwgLFxcL0JDU19wcm9wb3J0aW9uXzIwMjMucGRmXFwpLCBoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxuI3ByaW50KGJjcy5wcm9wLjIwMjMpXG4jZGV2Lm9mZigpXG5cbiNCQ1NfcHJvcG9ydGlvbl9hbGwrc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYyhcXCMwMDNDMzBcXCxcXCMwMTY2NUVcXCxcXCM4MENEQzFcXCxcXCNDN0VBRTVcXCxcXCNGNkU4QzNcXCxcXCNERkMyN0RcXCxcXCNCRjgxMkRcXCwgXFwjOEM1MTBBXFwsXFwjNTQzMDA1XFwgKSlcbiNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKFxcIzAwM0MzMFxcLFxcIzAxNjY1RVxcLFxcIzgwQ0RDMVxcLFxcI0M3RUFFNVxcLFxcYXp1cmUxXFwsXFxzbGF0ZWdyYXkxXFwsXFxzbGF0ZWdyYXkzXFwsIFxcc2xhdGVncmF5NFxcLCBcXGdyYXkyNVxcICkpIFxuXG5mIDwtIGZ1bmN0aW9uKHBhbCkgYnJld2VyLnBhbChicmV3ZXIucGFsLmluZm9bcGFsLCBcXG1heGNvbG9yc1xcXSwgcGFsKVxuKGNvbHMgPC0gZihcXFlsR25CdVxcKSlcbmBgYFxuXG48IS0tIHJuYi1zb3VyY2UtZW5kIC0tPlxuIn0= -->

#Coding graphs for LIS Conference - only 2023 data
data_2023 <- data_all%>%  filter(!year==\2024\)

df_BCS_proportions<- data_2023 %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

df_BCS_proportions<- na.omit(df_BCS_proportions)

df_BCS_proportions$month <- factor(df_BCS_proportions$month, levels = c(\5\,\6\,\7\, \8\, \9\, \10\, \11\),
        labels=c(\May\, \June\, \July\, \Aug\, \Sept\, \Oct\, \Nov\))
df_BCS_proportions$site <- factor(df_BCS_proportions$site, levels = c(\ASHC\,\FENC\,\GOLD\, \LAUR\),
        labels=c(\Ash Creek\, \Fence Creek\, \Gold Star\, \Laurel Hollow\))

bcs.prop.2023 <- ggplot(data=df_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat=\identity\, position = \fill\, colour = \black\)+  
   theme_bw() +  theme(panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1, size = rel(.9)))+
  theme(panel.grid.major.x =element_blank())+
  labs(title=\Body condition scores in 2023\, x =\Month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.2), angle =90), axis.title.x = element_text(size = rel(1.2), angle = 0))+ scale_fill_brewer(palette = \Blues\) +
  theme(strip.text = element_text(size = 15))+
  theme (legend.title = element_text(size = 18))+
  theme(title = element_text(size = 17))+
  theme(axis.text=element_text(size=17))+
  labs(fill =\Body Condition Score\)+
  facet_wrap(~ site )
bcs.prop.2023

#pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/BCS_proportion_2023.pdf\), height = 7, width = 13)
#print(bcs.prop.2023)
#dev.off()

#BCS_proportion_all+scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\#F6E8C3\,\#DFC27D\,\#BF812D\, \#8C510A\,\#543005\ ))
#scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\azure1\,\slategray1\,\slategray3\, \slategray4\, \gray25\ )) 

f <- function(pal) brewer.pal(brewer.pal.info[pal, \maxcolors\], pal)
(cols <- f(\YlGnBu\))

<!-- rnb-output-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dUkwTnZaR2x1WnlCbmNtRndhSE1nWm05eUlFeEpVeUJEYjI1bVpYSmxibU5sSUMwZ2IyNXNlU0F5TURJeklHUmhkR0ZjYm1SaGRHRmZNakF5TXlBOExTQmtZWFJoWDJGc2JDVStKU0FnWm1sc2RHVnlLQ0Y1WldGeVBUMWNYREl3TWpSY1hDbGNibHh1WkdaZlFrTlRYM0J5YjNCdmNuUnBiMjV6UEMwZ1pHRjBZVjh5TURJeklDVStKVnh1SUNCdGRYUmhkR1VvWTI5dVpHbDBhVzl1WDNOamIzSmxYMjUxYldWeWFXTWdQU0JoY3k1dWRXMWxjbWxqS0dOdmJtUnBkR2x2Ymw5elkyOXlaU2tzWTI5dVpHbDBhVzl1WDNOamIzSmxYMkpwYmlBOUlHTmhjMlZmZDJobGJpaGpiMjVrYVhScGIyNWZjMk52Y21VZ1BUMGdNU0IrSUZ4Y01WeGNMQ0JqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTWlCK1hGd3lYRndzWTI5dVpHbDBhVzl1WDNOamIzSmxJRDA5SURNZ2ZpQmNYRE5jWEN3Z1kyOXVaR2wwYVc5dVgzTmpiM0psSUQwOUlEUWdmbHhjTkZ4Y0xHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBTQTFJSDRnWEZ3MVhGd3NJR052Ym1ScGRHbHZibDl6WTI5eVpTQTlQU0EySUg1Y1hEWmNYQ3hqYjI1a2FYUnBiMjVmYzJOdmNtVWdQVDBnTnlCK0lGeGNOMXhjTENCamIyNWthWFJwYjI1ZmMyTnZjbVVnUFQwZ09DQitYRnc0WEZ3c0lHTnZibVJwZEdsdmJsOXpZMjl5WlNBOVBUa2dmbHhjT1Z4Y0xDQlVVbFZGSUg0Z1lYTXVZMmhoY21GamRHVnlLR052Ym1ScGRHbHZibDl6WTI5eVpTa3BLU0FsUGlWY2JpQWdaM0p2ZFhCZllua29jMmwwWlN3Z2JXOXVkR2dzSUdOdmJtUnBkR2x2Ymw5elkyOXlaVjlpYVc0cElDVStKU0JrY0d4NWNpQTZPbk4xYlcxaGNtbHpaU2hEYjNWdWREMGdiaWdwS1NBbFBpVmNiaUFnZFc1bmNtOTFjQ2dwSUNVK0pWeHVJQ0J0ZFhSaGRHVW9VSEp2Y0c5eWRHbHZiaUE5SUVOdmRXNTBMM04xYlNoRGIzVnVkQ2twWEc1Y2JtUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN6d3RJRzVoTG05dGFYUW9aR1pmUWtOVFgzQnliM0J2Y25ScGIyNXpLVnh1WEc1a1psOUNRMU5mY0hKdmNHOXlkR2x2Ym5Na2JXOXVkR2dnUEMwZ1ptRmpkRzl5S0dSbVgwSkRVMTl3Y205d2IzSjBhVzl1Y3lSdGIyNTBhQ3dnYkdWMlpXeHpJRDBnWXloY1hEVmNYQ3hjWERaY1hDeGNYRGRjWEN3Z1hGdzRYRndzSUZ4Y09WeGNMQ0JjWERFd1hGd3NJRnhjTVRGY1hDa3NYRzRnSUNBZ0lDQWdJR3hoWW1Wc2N6MWpLRnhjVFdGNVhGd3NJRnhjU25WdVpWeGNMQ0JjWEVwMWJIbGNYQ3dnWEZ4QmRXZGNYQ3dnWEZ4VFpYQjBYRndzSUZ4Y1QyTjBYRndzSUZ4Y1RtOTJYRndwS1Z4dVpHWmZRa05UWDNCeWIzQnZjblJwYjI1ekpITnBkR1VnUEMwZ1ptRmpkRzl5S0dSbVgwSkRVMTl3Y205d2IzSjBhVzl1Y3lSemFYUmxMQ0JzWlhabGJITWdQU0JqS0Z4Y1FWTklRMXhjTEZ4Y1JrVk9RMXhjTEZ4Y1IwOU1SRnhjTENCY1hFeEJWVkpjWENrc1hHNGdJQ0FnSUNBZ0lHeGhZbVZzY3oxaktGeGNRWE5vSUVOeVpXVnJYRndzSUZ4Y1JtVnVZMlVnUTNKbFpXdGNYQ3dnWEZ4SGIyeGtJRk4wWVhKY1hDd2dYRnhNWVhWeVpXd2dTRzlzYkc5M1hGd3BLVnh1WEc1aVkzTXVjSEp2Y0M0eU1ESXpJRHd0SUdkbmNHeHZkQ2hrWVhSaFBXUm1YMEpEVTE5d2NtOXdiM0owYVc5dWN5d2dZV1Z6S0hnOWJXOXVkR2dzSUhrOUlGQnliM0J2Y25ScGIyNHNJR1pwYkd3OVkyOXVaR2wwYVc5dVgzTmpiM0psWDJKcGJpa3BJQ3RjYmlBZ1oyVnZiVjlpWVhJb2QybGtkR2dnUFNBdU5Td2djM1JoZEQxY1hHbGtaVzUwYVhSNVhGd3NJSEJ2YzJsMGFXOXVJRDBnWEZ4bWFXeHNYRndzSUdOdmJHOTFjaUE5SUZ4Y1lteGhZMnRjWENrcklDQmNiaUFnSUhSb1pXMWxYMkozS0NrZ0t5QWdkR2hsYldVb2NHRnVaV3d1WjNKcFpDNXRhVzV2Y2lBOUlHVnNaVzFsYm5SZllteGhibXNvS1NrcklGeHVJQ0IwYUdWdFpTaGhlR2x6TG5SbGVIUXVlQ0E5SUdWc1pXMWxiblJmZEdWNGRDaGhibWRzWlNBOUlEUTFMQ0IyYW5WemRDQTlJREVzSUdocWRYTjBQVEVzSUhOcGVtVWdQU0J5Wld3b0xqa3BLU2tyWEc0Z0lIUm9aVzFsS0hCaGJtVnNMbWR5YVdRdWJXRnFiM0l1ZUNBOVpXeGxiV1Z1ZEY5aWJHRnVheWdwS1N0Y2JpQWdiR0ZpY3loMGFYUnNaVDFjWEVKdlpIa2dZMjl1WkdsMGFXOXVJSE5qYjNKbGN5QnBiaUF5TURJelhGd3NJSGdnUFZ4Y1RXOXVkR2hjWEN3Z2VTQTlJRnhjVUhKdmNHOXlkR2x2YmlCdlppQmliMlI1SUdOdmJtUnBkR2x2YmlCelkyOXlaWE1nWEZ3cEt5QjBhR1Z0WlNoaGVHbHpMblJwZEd4bExua2dQU0JsYkdWdFpXNTBYM1JsZUhRb2MybDZaU0E5SUhKbGJDZ3hMaklwTENCaGJtZHNaU0E5T1RBcExDQmhlR2x6TG5ScGRHeGxMbmdnUFNCbGJHVnRaVzUwWDNSbGVIUW9jMmw2WlNBOUlISmxiQ2d4TGpJcExDQmhibWRzWlNBOUlEQXBLU3NnYzJOaGJHVmZabWxzYkY5aWNtVjNaWElvY0dGc1pYUjBaU0E5SUZ4Y1FteDFaWE5jWENrZ0sxeHVJQ0IwYUdWdFpTaHpkSEpwY0M1MFpYaDBJRDBnWld4bGJXVnVkRjkwWlhoMEtITnBlbVVnUFNBeE5Ta3BLMXh1SUNCMGFHVnRaU0FvYkdWblpXNWtMblJwZEd4bElEMGdaV3hsYldWdWRGOTBaWGgwS0hOcGVtVWdQU0F4T0NrcEsxeHVJQ0IwYUdWdFpTaDBhWFJzWlNBOUlHVnNaVzFsYm5SZmRHVjRkQ2h6YVhwbElEMGdNVGNwS1N0Y2JpQWdkR2hsYldVb1lYaHBjeTUwWlhoMFBXVnNaVzFsYm5SZmRHVjRkQ2h6YVhwbFBURTNLU2tyWEc0Z0lHeGhZbk1vWm1sc2JDQTlYRnhDYjJSNUlFTnZibVJwZEdsdmJpQlRZMjl5WlZ4Y0tTdGNiaUFnWm1GalpYUmZkM0poY0NoK0lITnBkR1VnS1Z4dVltTnpMbkJ5YjNBdU1qQXlNMXh1WEc0amNHUm1LSEJoYzNSbE1DaHdZWFJvSUQwZ1hGeE1ZV0pmUkdGMFlWOVVhWE56ZFdWUWNtOWpaWE56YVc1bkwyOTFkSEIxZEZ4Y0lDeGNYQzlDUTFOZmNISnZjRzl5ZEdsdmJsOHlNREl6TG5Ca1pseGNLU3dnYUdWcFoyaDBJRDBnTnl3Z2QybGtkR2dnUFNBeE15bGNiaU53Y21sdWRDaGlZM011Y0hKdmNDNHlNREl6S1Z4dUkyUmxkaTV2Wm1Zb0tWeHVYRzRqUWtOVFgzQnliM0J2Y25ScGIyNWZZV3hzSzNOallXeGxYMlpwYkd4ZmJXRnVkV0ZzS0haaGJIVmxjeUE5SUdNb1hGd2pNREF6UXpNd1hGd3NYRndqTURFMk5qVkZYRndzWEZ3ak9EQkRSRU14WEZ3c1hGd2pRemRGUVVVMVhGd3NYRndqUmpaRk9FTXpYRndzWEZ3alJFWkRNamRFWEZ3c1hGd2pRa1k0TVRKRVhGd3NJRnhjSXpoRE5URXdRVnhjTEZ4Y0l6VTBNekF3TlZ4Y0lDa3BYRzRqYzJOaGJHVmZabWxzYkY5dFlXNTFZV3dvZG1Gc2RXVnpJRDBnWXloY1hDTXdNRE5ETXpCY1hDeGNYQ013TVRZMk5VVmNYQ3hjWENNNE1FTkVRekZjWEN4Y1hDTkROMFZCUlRWY1hDeGNYR0Y2ZFhKbE1WeGNMRnhjYzJ4aGRHVm5jbUY1TVZ4Y0xGeGNjMnhoZEdWbmNtRjVNMXhjTENCY1hITnNZWFJsWjNKaGVUUmNYQ3dnWEZ4bmNtRjVNalZjWENBcEtTQmNibHh1WmlBOExTQm1kVzVqZEdsdmJpaHdZV3dwSUdKeVpYZGxjaTV3WVd3b1luSmxkMlZ5TG5CaGJDNXBibVp2VzNCaGJDd2dYRnh0WVhoamIyeHZjbk5jWEYwc0lIQmhiQ2xjYmloamIyeHpJRHd0SUdZb1hGeFpiRWR1UW5WY1hDa3BYRzVnWUdCY2JtQmdZQ0o5IC0tPlxuXG5gYGByXG5gYGByXG4jQ29kaW5nIGdyYXBocyBmb3IgTElTIENvbmZlcmVuY2UgLSBvbmx5IDIwMjMgZGF0YVxuZGF0YV8yMDIzIDwtIGRhdGFfYWxsJT4lICBmaWx0ZXIoIXllYXI9PVxcMjAyNFxcKVxuXG5kZl9CQ1NfcHJvcG9ydGlvbnM8LSBkYXRhXzIwMjMgJT4lXG4gIG11dGF0ZShjb25kaXRpb25fc2NvcmVfbnVtZXJpYyA9IGFzLm51bWVyaWMoY29uZGl0aW9uX3Njb3JlKSxjb25kaXRpb25fc2NvcmVfYmluID0gY2FzZV93aGVuKGNvbmRpdGlvbl9zY29yZSA9PSAxIH4gXFwxXFwsIGNvbmRpdGlvbl9zY29yZSA9PSAyIH5cXDJcXCxjb25kaXRpb25fc2NvcmUgPT0gMyB+IFxcM1xcLCBjb25kaXRpb25fc2NvcmUgPT0gNCB+XFw0XFwsY29uZGl0aW9uX3Njb3JlID09IDUgfiBcXDVcXCwgY29uZGl0aW9uX3Njb3JlID09IDYgflxcNlxcLGNvbmRpdGlvbl9zY29yZSA9PSA3IH4gXFw3XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA4IH5cXDhcXCwgY29uZGl0aW9uX3Njb3JlID09OSB+XFw5XFwsIFRSVUUgfiBhcy5jaGFyYWN0ZXIoY29uZGl0aW9uX3Njb3JlKSkpICU+JVxuICBncm91cF9ieShzaXRlLCBtb250aCwgY29uZGl0aW9uX3Njb3JlX2JpbikgJT4lIGRwbHlyIDo6c3VtbWFyaXNlKENvdW50PSBuKCkpICU+JVxuICB1bmdyb3VwKCkgJT4lXG4gIG11dGF0ZShQcm9wb3J0aW9uID0gQ291bnQvc3VtKENvdW50KSlcblxuZGZfQkNTX3Byb3BvcnRpb25zPC0gbmEub21pdChkZl9CQ1NfcHJvcG9ydGlvbnMpXG5cbmRmX0JDU19wcm9wb3J0aW9ucyRtb250aCA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJG1vbnRoLCBsZXZlbHMgPSBjKFxcNVxcLFxcNlxcLFxcN1xcLCBcXDhcXCwgXFw5XFwsIFxcMTBcXCwgXFwxMVxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxNYXlcXCwgXFxKdW5lXFwsIFxcSnVseVxcLCBcXEF1Z1xcLCBcXFNlcHRcXCwgXFxPY3RcXCwgXFxOb3ZcXCkpXG5kZl9CQ1NfcHJvcG9ydGlvbnMkc2l0ZSA8LSBmYWN0b3IoZGZfQkNTX3Byb3BvcnRpb25zJHNpdGUsIGxldmVscyA9IGMoXFxBU0hDXFwsXFxGRU5DXFwsXFxHT0xEXFwsIFxcTEFVUlxcKSxcbiAgICAgICAgbGFiZWxzPWMoXFxBc2ggQ3JlZWtcXCwgXFxGZW5jZSBDcmVla1xcLCBcXEdvbGQgU3RhclxcLCBcXExhdXJlbCBIb2xsb3dcXCkpXG5cbmJjcy5wcm9wLjIwMjMgPC0gZ2dwbG90KGRhdGE9ZGZfQkNTX3Byb3BvcnRpb25zLCBhZXMoeD1tb250aCwgeT0gUHJvcG9ydGlvbiwgZmlsbD1jb25kaXRpb25fc2NvcmVfYmluKSkgK1xuICBnZW9tX2Jhcih3aWR0aCA9IC41LCBzdGF0PVxcaWRlbnRpdHlcXCwgcG9zaXRpb24gPSBcXGZpbGxcXCwgY29sb3VyID0gXFxibGFja1xcKSsgIFxuICAgdGhlbWVfYncoKSArICB0aGVtZShwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpKSsgXG4gIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIHZqdXN0ID0gMSwgaGp1c3Q9MSwgc2l6ZSA9IHJlbCguOSkpKStcbiAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvci54ID1lbGVtZW50X2JsYW5rKCkpK1xuICBsYWJzKHRpdGxlPVxcQm9keSBjb25kaXRpb24gc2NvcmVzIGluIDIwMjNcXCwgeCA9XFxNb250aFxcLCB5ID0gXFxQcm9wb3J0aW9uIG9mIGJvZHkgY29uZGl0aW9uIHNjb3JlcyBcXCkrIHRoZW1lKGF4aXMudGl0bGUueSA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID05MCksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gcmVsKDEuMiksIGFuZ2xlID0gMCkpKyBzY2FsZV9maWxsX2JyZXdlcihwYWxldHRlID0gXFxCbHVlc1xcKSArXG4gIHRoZW1lKHN0cmlwLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkrXG4gIHRoZW1lIChsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE4KSkrXG4gIHRoZW1lKHRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNykpK1xuICB0aGVtZShheGlzLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9MTcpKStcbiAgbGFicyhmaWxsID1cXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwpK1xuICBmYWNldF93cmFwKH4gc2l0ZSApXG5iY3MucHJvcC4yMDIzXG5cbiNwZGYocGFzdGUwKHBhdGggPSBcXExhYl9EYXRhX1Rpc3N1ZVByb2Nlc3Npbmcvb3V0cHV0XFwgLFxcL0JDU19wcm9wb3J0aW9uXzIwMjMucGRmXFwpLCBoZWlnaHQgPSA3LCB3aWR0aCA9IDEzKVxuI3ByaW50KGJjcy5wcm9wLjIwMjMpXG4jZGV2Lm9mZigpXG5cbiNCQ1NfcHJvcG9ydGlvbl9hbGwrc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYyhcXCMwMDNDMzBcXCxcXCMwMTY2NUVcXCxcXCM4MENEQzFcXCxcXCNDN0VBRTVcXCxcXCNGNkU4QzNcXCxcXCNERkMyN0RcXCxcXCNCRjgxMkRcXCwgXFwjOEM1MTBBXFwsXFwjNTQzMDA1XFwgKSlcbiNzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKFxcIzAwM0MzMFxcLFxcIzAxNjY1RVxcLFxcIzgwQ0RDMVxcLFxcI0M3RUFFNVxcLFxcYXp1cmUxXFwsXFxzbGF0ZWdyYXkxXFwsXFxzbGF0ZWdyYXkzXFwsIFxcc2xhdGVncmF5NFxcLCBcXGdyYXkyNVxcICkpIFxuXG5mIDwtIGZ1bmN0aW9uKHBhbCkgYnJld2VyLnBhbChicmV3ZXIucGFsLmluZm9bcGFsLCBcXG1heGNvbG9yc1xcXSwgcGFsKVxuKGNvbHMgPC0gZihcXFlsR25CdVxcKSlcbmBgYFxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuI0NvZGluZyBncmFwaHMgZm9yIExJUyBDb25mZXJlbmNlIC0gb25seSAyMDIzIGRhdGFcbmRhdGFfMjAyMyA8LSBkYXRhX2FsbCU+JSAgZmlsdGVyKCF5ZWFyPT1cXDIwMjRcXClcblxuZGZfQkNTX3Byb3BvcnRpb25zPC0gZGF0YV8yMDIzICU+JVxuICBtdXRhdGUoY29uZGl0aW9uX3Njb3JlX251bWVyaWMgPSBhcy5udW1lcmljKGNvbmRpdGlvbl9zY29yZSksY29uZGl0aW9uX3Njb3JlX2JpbiA9IGNhc2Vfd2hlbihjb25kaXRpb25fc2NvcmUgPT0gMSB+IFxcMVxcLCBjb25kaXRpb25fc2NvcmUgPT0gMiB+XFwyXFwsY29uZGl0aW9uX3Njb3JlID09IDMgfiBcXDNcXCwgY29uZGl0aW9uX3Njb3JlID09IDQgflxcNFxcLGNvbmRpdGlvbl9zY29yZSA9PSA1IH4gXFw1XFwsIGNvbmRpdGlvbl9zY29yZSA9PSA2IH5cXDZcXCxjb25kaXRpb25fc2NvcmUgPT0gNyB+IFxcN1xcLCBjb25kaXRpb25fc2NvcmUgPT0gOCB+XFw4XFwsIGNvbmRpdGlvbl9zY29yZSA9PTkgflxcOVxcLCBUUlVFIH4gYXMuY2hhcmFjdGVyKGNvbmRpdGlvbl9zY29yZSkpKSAlPiVcbiAgZ3JvdXBfYnkoc2l0ZSwgbW9udGgsIGNvbmRpdGlvbl9zY29yZV9iaW4pICU+JSBkcGx5ciA6OnN1bW1hcmlzZShDb3VudD0gbigpKSAlPiVcbiAgdW5ncm91cCgpICU+JVxuICBtdXRhdGUoUHJvcG9ydGlvbiA9IENvdW50L3N1bShDb3VudCkpXG5cbmRmX0JDU19wcm9wb3J0aW9uczwtIG5hLm9taXQoZGZfQkNTX3Byb3BvcnRpb25zKVxuXG5kZl9CQ1NfcHJvcG9ydGlvbnMkbW9udGggPC0gZmFjdG9yKGRmX0JDU19wcm9wb3J0aW9ucyRtb250aCwgbGV2ZWxzID0gYyhcXDVcXCxcXDZcXCxcXDdcXCwgXFw4XFwsIFxcOVxcLCBcXDEwXFwsIFxcMTFcXCksXG4gICAgICAgIGxhYmVscz1jKFxcTWF5XFwsIFxcSnVuZVxcLCBcXEp1bHlcXCwgXFxBdWdcXCwgXFxTZXB0XFwsIFxcT2N0XFwsIFxcTm92XFwpKVxuZGZfQkNTX3Byb3BvcnRpb25zJHNpdGUgPC0gZmFjdG9yKGRmX0JDU19wcm9wb3J0aW9ucyRzaXRlLCBsZXZlbHMgPSBjKFxcQVNIQ1xcLFxcRkVOQ1xcLFxcR09MRFxcLCBcXExBVVJcXCksXG4gICAgICAgIGxhYmVscz1jKFxcQXNoIENyZWVrXFwsIFxcRmVuY2UgQ3JlZWtcXCwgXFxHb2xkIFN0YXJcXCwgXFxMYXVyZWwgSG9sbG93XFwpKVxuXG5iY3MucHJvcC4yMDIzIDwtIGdncGxvdChkYXRhPWRmX0JDU19wcm9wb3J0aW9ucywgYWVzKHg9bW9udGgsIHk9IFByb3BvcnRpb24sIGZpbGw9Y29uZGl0aW9uX3Njb3JlX2JpbikpICtcbiAgZ2VvbV9iYXIod2lkdGggPSAuNSwgc3RhdD1cXGlkZW50aXR5XFwsIHBvc2l0aW9uID0gXFxmaWxsXFwsIGNvbG91ciA9IFxcYmxhY2tcXCkrICBcbiAgIHRoZW1lX2J3KCkgKyAgdGhlbWUocGFuZWwuZ3JpZC5taW5vciA9IGVsZW1lbnRfYmxhbmsoKSkrIFxuICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDQ1LCB2anVzdCA9IDEsIGhqdXN0PTEsIHNpemUgPSByZWwoLjkpKSkrXG4gIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IueCA9ZWxlbWVudF9ibGFuaygpKStcbiAgbGFicyh0aXRsZT1cXEJvZHkgY29uZGl0aW9uIHNjb3JlcyBpbiAyMDIzXFwsIHggPVxcTW9udGhcXCwgeSA9IFxcUHJvcG9ydGlvbiBvZiBib2R5IGNvbmRpdGlvbiBzY29yZXMgXFwpKyB0aGVtZShheGlzLnRpdGxlLnkgPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjIpLCBhbmdsZSA9OTApLCBheGlzLnRpdGxlLnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IHJlbCgxLjIpLCBhbmdsZSA9IDApKSsgc2NhbGVfZmlsbF9icmV3ZXIocGFsZXR0ZSA9IFxcQmx1ZXNcXCkgK1xuICB0aGVtZShzdHJpcC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSkpK1xuICB0aGVtZSAobGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxOCkpK1xuICB0aGVtZSh0aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTcpKStcbiAgdGhlbWUoYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTE3KSkrXG4gIGxhYnMoZmlsbCA9XFxCb2R5IENvbmRpdGlvbiBTY29yZVxcKStcbiAgZmFjZXRfd3JhcCh+IHNpdGUgKVxuYmNzLnByb3AuMjAyM1xuXG4jcGRmKHBhc3RlMChwYXRoID0gXFxMYWJfRGF0YV9UaXNzdWVQcm9jZXNzaW5nL291dHB1dFxcICxcXC9CQ1NfcHJvcG9ydGlvbl8yMDIzLnBkZlxcKSwgaGVpZ2h0ID0gNywgd2lkdGggPSAxMylcbiNwcmludChiY3MucHJvcC4yMDIzKVxuI2Rldi5vZmYoKVxuXG4jQkNTX3Byb3BvcnRpb25fYWxsK3NjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoXFwjMDAzQzMwXFwsXFwjMDE2NjVFXFwsXFwjODBDREMxXFwsXFwjQzdFQUU1XFwsXFwjRjZFOEMzXFwsXFwjREZDMjdEXFwsXFwjQkY4MTJEXFwsIFxcIzhDNTEwQVxcLFxcIzU0MzAwNVxcICkpXG4jc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYyhcXCMwMDNDMzBcXCxcXCMwMTY2NUVcXCxcXCM4MENEQzFcXCxcXCNDN0VBRTVcXCxcXGF6dXJlMVxcLFxcc2xhdGVncmF5MVxcLFxcc2xhdGVncmF5M1xcLCBcXHNsYXRlZ3JheTRcXCwgXFxncmF5MjVcXCApKSBcblxuZiA8LSBmdW5jdGlvbihwYWwpIGJyZXdlci5wYWwoYnJld2VyLnBhbC5pbmZvW3BhbCwgXFxtYXhjb2xvcnNcXF0sIHBhbClcbihjb2xzIDwtIGYoXFxZbEduQnVcXCkpXG5gYGBcbmBgYCJ9 -->

```r
```r
#Coding graphs for LIS Conference - only 2023 data
data_2023 <- data_all%>%  filter(!year==\2024\)

df_BCS_proportions<- data_2023 %>%
  mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ \1\, condition_score == 2 ~\2\,condition_score == 3 ~ \3\, condition_score == 4 ~\4\,condition_score == 5 ~ \5\, condition_score == 6 ~\6\,condition_score == 7 ~ \7\, condition_score == 8 ~\8\, condition_score ==9 ~\9\, TRUE ~ as.character(condition_score))) %>%
  group_by(site, month, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
  ungroup() %>%
  mutate(Proportion = Count/sum(Count))

df_BCS_proportions<- na.omit(df_BCS_proportions)

df_BCS_proportions$month <- factor(df_BCS_proportions$month, levels = c(\5\,\6\,\7\, \8\, \9\, \10\, \11\),
        labels=c(\May\, \June\, \July\, \Aug\, \Sept\, \Oct\, \Nov\))
df_BCS_proportions$site <- factor(df_BCS_proportions$site, levels = c(\ASHC\,\FENC\,\GOLD\, \LAUR\),
        labels=c(\Ash Creek\, \Fence Creek\, \Gold Star\, \Laurel Hollow\))

bcs.prop.2023 <- ggplot(data=df_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
  geom_bar(width = .5, stat=\identity\, position = \fill\, colour = \black\)+  
   theme_bw() +  theme(panel.grid.minor = element_blank())+ 
  theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1, size = rel(.9)))+
  theme(panel.grid.major.x =element_blank())+
  labs(title=\Body condition scores in 2023\, x =\Month\, y = \Proportion of body condition scores \)+ theme(axis.title.y = element_text(size = rel(1.2), angle =90), axis.title.x = element_text(size = rel(1.2), angle = 0))+ scale_fill_brewer(palette = \Blues\) +
  theme(strip.text = element_text(size = 15))+
  theme (legend.title = element_text(size = 18))+
  theme(title = element_text(size = 17))+
  theme(axis.text=element_text(size=17))+
  labs(fill =\Body Condition Score\)+
  facet_wrap(~ site )
bcs.prop.2023

#pdf(paste0(path = \Lab_Data_TissueProcessing/output\ ,\/BCS_proportion_2023.pdf\), height = 7, width = 13)
#print(bcs.prop.2023)
#dev.off()

#BCS_proportion_all+scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\#F6E8C3\,\#DFC27D\,\#BF812D\, \#8C510A\,\#543005\ ))
#scale_fill_manual(values = c(\#003C30\,\#01665E\,\#80CDC1\,\#C7EAE5\,\azure1\,\slategray1\,\slategray3\, \slategray4\, \gray25\ )) 

f <- function(pal) brewer.pal(brewer.pal.info[pal, \maxcolors\], pal)
(cols <- f(\YlGnBu\))
```
```

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-text-begin -->



<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVJME52WkdsdVp5Qm5jbUZ3YUhNZ1ptOXlJRXhKVXlCRGIyNW1aWEpsYm1ObElDMGdiMjVzZVNBeU1ESXpJR1JoZEdFZ0tHTnZiblJwYm5WbFpDbGNibHh1WkdaZk0xOXdjbTl3UEMwZ2MyVnNaV04wS0dSaGRHRmZNakF5TXl3Z0xXTW9YQ0p6YVhSbFhDSXBLVnh1WEc0alEyOTFiblFnWVc1a0lIQnliM0J2Y25ScGIyNGdkR0ZpYkdVZ1ptOXlJREl3TWpNZ1ltTnpJSE5qYjNKbGN5QmllU0J0YjI1MGFGeHVZbU56WDNSaFlteGxJRHd0SUhSaFlteGxLR1JoZEdGZk1qQXlNeVJqYjI1a2FYUnBiMjVmYzJOdmNtVXNJR1JoZEdGZk1qQXlNeVJ0YjI1MGFDd2daR0YwWVY4eU1ESXpKSE5wZEdVcFhHNXVZVzFsY3loa2FXMXVZVzFsY3loaVkzTmZkR0ZpYkdVcEtTQThMU0JqS0Z3aVFtOWtlU0JEYjI1a2FYUnBiMjRnVTJOdmNtVmNJaXdnWENKTmIyNTBhRndpTENCY0lsTnBkR1ZjSWlsY2JtSmpjMTkwWVdKc1pTQThMU0JoWkdSdFlYSm5hVzV6S0dKamMxOTBZV0pzWlNsY2JtSmpjMTkwWVdKc1pWeHVYRzVnWUdBaWZRPT0gLS0+XG5cbmBgYHJcbiNDb2RpbmcgZ3JhcGhzIGZvciBMSVMgQ29uZmVyZW5jZSAtIG9ubHkgMjAyMyBkYXRhIChjb250aW51ZWQpXG5cbmRmXzNfcHJvcDwtIHNlbGVjdChkYXRhXzIwMjMsIC1jKFxcc2l0ZVxcKSlcblxuI0NvdW50IGFuZCBwcm9wb3J0aW9uIHRhYmxlIGZvciAyMDIzIGJjcyBzY29yZXMgYnkgbW9udGhcbmJjc190YWJsZSA8LSB0YWJsZShkYXRhXzIwMjMkY29uZGl0aW9uX3Njb3JlLCBkYXRhXzIwMjMkbW9udGgsIGRhdGFfMjAyMyRzaXRlKVxubmFtZXMoZGltbmFtZXMoYmNzX3RhYmxlKSkgPC0gYyhcXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwsIFxcTW9udGhcXCwgXFxTaXRlXFwpXG5iY3NfdGFibGUgPC0gYWRkbWFyZ2lucyhiY3NfdGFibGUpXG5iY3NfdGFibGVcblxuYGBgXG5cbjwhLS0gcm5iLXNvdXJjZS1lbmQgLS0+XG4ifQ== -->

#Coding graphs for LIS Conference - only 2023 data (continued)

df_3_prop<- select(data_2023, -c(\site\))

#Count and proportion table for 2023 bcs scores by month
bcs_table <- table(data_2023$condition_score, data_2023$month, data_2023$site)
names(dimnames(bcs_table)) <- c(\Body Condition Score\, \Month\, \Site\)
bcs_table <- addmargins(bcs_table)
bcs_table

<!-- rnb-output-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVZR0JnY2x4dUkwTnZaR2x1WnlCbmNtRndhSE1nWm05eUlFeEpVeUJEYjI1bVpYSmxibU5sSUMwZ2IyNXNlU0F5TURJeklHUmhkR0VnS0dOdmJuUnBiblZsWkNsY2JseHVaR1pmTTE5d2NtOXdQQzBnYzJWc1pXTjBLR1JoZEdGZk1qQXlNeXdnTFdNb1hGeHphWFJsWEZ3cEtWeHVYRzRqUTI5MWJuUWdZVzVrSUhCeWIzQnZjblJwYjI0Z2RHRmliR1VnWm05eUlESXdNak1nWW1OeklITmpiM0psY3lCaWVTQnRiMjUwYUZ4dVltTnpYM1JoWW14bElEd3RJSFJoWW14bEtHUmhkR0ZmTWpBeU15UmpiMjVrYVhScGIyNWZjMk52Y21Vc0lHUmhkR0ZmTWpBeU15UnRiMjUwYUN3Z1pHRjBZVjh5TURJekpITnBkR1VwWEc1dVlXMWxjeWhrYVcxdVlXMWxjeWhpWTNOZmRHRmliR1VwS1NBOExTQmpLRnhjUW05a2VTQkRiMjVrYVhScGIyNGdVMk52Y21WY1hDd2dYRnhOYjI1MGFGeGNMQ0JjWEZOcGRHVmNYQ2xjYm1KamMxOTBZV0pzWlNBOExTQmhaR1J0WVhKbmFXNXpLR0pqYzE5MFlXSnNaU2xjYm1KamMxOTBZV0pzWlZ4dVhHNWdZR0JjYm1CZ1lDSjkgLS0+XG5cbmBgYHJcbmBgYHJcbiNDb2RpbmcgZ3JhcGhzIGZvciBMSVMgQ29uZmVyZW5jZSAtIG9ubHkgMjAyMyBkYXRhIChjb250aW51ZWQpXG5cbmRmXzNfcHJvcDwtIHNlbGVjdChkYXRhXzIwMjMsIC1jKFxcc2l0ZVxcKSlcblxuI0NvdW50IGFuZCBwcm9wb3J0aW9uIHRhYmxlIGZvciAyMDIzIGJjcyBzY29yZXMgYnkgbW9udGhcbmJjc190YWJsZSA8LSB0YWJsZShkYXRhXzIwMjMkY29uZGl0aW9uX3Njb3JlLCBkYXRhXzIwMjMkbW9udGgsIGRhdGFfMjAyMyRzaXRlKVxubmFtZXMoZGltbmFtZXMoYmNzX3RhYmxlKSkgPC0gYyhcXEJvZHkgQ29uZGl0aW9uIFNjb3JlXFwsIFxcTW9udGhcXCwgXFxTaXRlXFwpXG5iY3NfdGFibGUgPC0gYWRkbWFyZ2lucyhiY3NfdGFibGUpXG5iY3NfdGFibGVcblxuYGBgXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuYGBgclxuI0NvZGluZyBncmFwaHMgZm9yIExJUyBDb25mZXJlbmNlIC0gb25seSAyMDIzIGRhdGEgKGNvbnRpbnVlZClcblxuZGZfM19wcm9wPC0gc2VsZWN0KGRhdGFfMjAyMywgLWMoXFxzaXRlXFwpKVxuXG4jQ291bnQgYW5kIHByb3BvcnRpb24gdGFibGUgZm9yIDIwMjMgYmNzIHNjb3JlcyBieSBtb250aFxuYmNzX3RhYmxlIDwtIHRhYmxlKGRhdGFfMjAyMyRjb25kaXRpb25fc2NvcmUsIGRhdGFfMjAyMyRtb250aCwgZGF0YV8yMDIzJHNpdGUpXG5uYW1lcyhkaW1uYW1lcyhiY3NfdGFibGUpKSA8LSBjKFxcQm9keSBDb25kaXRpb24gU2NvcmVcXCwgXFxNb250aFxcLCBcXFNpdGVcXClcbmJjc190YWJsZSA8LSBhZGRtYXJnaW5zKGJjc190YWJsZSlcbmJjc190YWJsZVxuXG5gYGBcbmBgYCJ9 -->

```r
```r
#Coding graphs for LIS Conference - only 2023 data (continued)

df_3_prop<- select(data_2023, -c(\site\))

#Count and proportion table for 2023 bcs scores by month
bcs_table <- table(data_2023$condition_score, data_2023$month, data_2023$site)
names(dimnames(bcs_table)) <- c(\Body Condition Score\, \Month\, \Site\)
bcs_table <- addmargins(bcs_table)
bcs_table

```
```

<!-- rnb-source-end -->


<!-- rnb-output-end -->

<!-- rnb-chunk-end -->


<!-- rnb-text-begin -->




<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHVjMlZ6YzJsdmJrbHVabThvS1Z4dVlHQmdJbjA9IC0tPlxuXG5gYGByXG5zZXNzaW9uSW5mbygpXG5gYGBcblxuPCEtLSBybmItc291cmNlLWVuZCAtLT5cbiJ9 -->

sessionInfo()

````

```r
sessionInfo()

```

Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.

When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).

The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike Knit, Preview does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.